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SUMMARY 


^ An  efficient  algorithm  for  a vortex  panel  method  in  the  complex 
plane  is  developed  to  compute  the  potential  flow  past  multi-component 
airfoils  in  free  air  and  a porous-wall  wind  tunnel.  The  theoretical  foundations 
of  the  method  — the  properties  of  the  source  and  vortex  density  functions, 
the  relationship  between  the  exterior  and  interior  flows,  and  the  general 
Kutta-Joukowski  condition  for  a trailing  edge  with  crossflow  — are  derived 
from  the  theory  of  the  Cauchy  type  integral.  The  method  utilizes  flat  panels 
with  linear  vortex  and  source  densities,  the  latter  being  used  to  simulate  the 
displacement  effect  of  boundary  layers.  The  airfoil  boundary  condition  is 
satisfied  at  all  panel  midpoints  and  the  overdetermined  system  of  linear 
algebraic  equations  solved  as  a least  squares  problem,  in  the(£j\norm,  or  in 
the  L^snorm.  The  wind  tunnel  wall  interference  problem  is  treated  using  the 
concept  of  the  Green’s  function  in  the  complex  plane.  Examples  are  worked 
out  for  some  theoretical  airfoils  and  extensive  tables  of  exact  and  computed 
pressure  distributions  are  given. 


L s u L 


RESUME 


Le  present  ouvrage  porte  sur  la  mise  au  point  d’un  algorithme 
efficace  qui  permettrait  de  calculer  dans  le  plan  complexe,  a partir  de  la  me- 
thode  des  petits  paves  a tourbillons,  l’ecoulement  a potentiel  autour  des 
profils  a composants  multiples  dans  l’air  libre  et  dans  une  soufflerie  a paroi 
permeable.  Les  fondements  theoriques  de  la  methode  — les  proprietes  des 
fonctions  de  densite  des  tourbillons  et  des  sources,  les  relations  entre  l’ecou- 
lement  exterieur  et  interieur  et  la  condition  generate  de  Kutta-Joukowski 
pour  un  bord  de  fuite  avec  ecoulement  transversal  — sont  tires  de  la  theorie 
de  l’integrale  de  Cauchy.  La  methode  fait  intervenir  des  paves  plans  presen- 
tant  des  fonctions  lineaires  de  densite  des  tourbillons  et  des  sources,  ces 
demieres  permettant  de  simuler  l’effet  de  deplacement  des  couches  limites. 
Les  conditions  aux  limites  des  profils  sont  satisfaites  en  tout  point  median 
des  paves,  et  le  systeme  algebrique  d’equations  lineaires  hyperdeterminees 
est  resolu  par  la  methode  des  moindres  carres  ou  de  la  norme  L)  ou  L„ . 
L’effet  de  la  paroi  de  la  soufflerie  est  traite  en  termes  du  concept  de  la  fonc- 


tion  de  Green  dans  le  plan  complexe.  Quelques  exemples  de  profils  theoriques 
sont  analyses  et  des  tables  completes  des  valeurs  exactes  de  differentes  re- 
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CALCULATION  OF  THE  POTENTIAL  FLOW  PAST  MULTI  COMPONENT  AIRFOILS 
USING  A VORTEX  PANEL  METHOD  IN  THE  COMPLEX  PLANE 


1.0  INTRODUCTION 

The  simplicity  and  compactness  of  a two-dimensional  panel  method  developed  in  the  complex 
plane  was  first  demonstrated  by  Giesing,  Reference  1,  on  flat  panels  with  a constant  source  density, 
which  are  the  building  stones  of  the  Douglas-Neumann  source  distribution  method.  Having  in  mind  a 
similar  approach  to  a contemporary  vortex  panel  method,  which  uses  linear  distributions  of  source  and 
vortex  densities  over  the  panels,  the  present  report  outlines  first  a theory  of  distributions  of  sources  and 
vortices  in  the  complex  plane.  The  properties  of  the  source  and  vortex  density  functions  and  the  rela- 
tionship between  the  exterior  and  interior  flows  are  derived  using  the  concept  of  the  Cauchy  type 
integral.  The  novelty  of  the  present  approach  is  that  the  integral  is  linked  to  the  complex  disturbance 
velocity  and  not  to  the  potential  function,  so  that  no  artificial  discontinuity  lines  need  to  be  introduced 
in  the  flow  field.  The  Cauchy  type  integral  is  also  used  to  derive  the  general  Kutta-Joukowski  condition 
in  terms  of  the  trailing  edge  values  of  the  source  and  vortex  densities.  The  connection  between  the 
conditions  of  a finite  trailing  edge  velocity  and  the  geometrical  concept  of  equal  upper  and  lower 
trailing  edge  velocities  is  elucidated  upon. 

The  major  part  of  the  work  is  focussed  on  the  development  of  an  algorithm  for  a vortex  panel 
method  which  uses  flat  panels  with  linear  distributions  of  source  and  vortex  densities.  The  comer  values 
of  the  source  density,  used  to  simulate  the  displacement  effect  of  boundary  layers,  enter  as  given  quan- 
tities and  the  comer  values  of  the  vortex  density,  determining  tangent  components  of  velocity,  are  the 
quantities  to  be  solved  for.  Since  the  airfoil  boundary  condition  is  met  at  all  panel  midpoints  and  the 
Kutta-Joukowski  condition  is  satisfied  at  trailing  edges,  the  resultant  system  of  linear  algebraic  equa- 
tions is  overdetermined  and  solved  by  minimizing  the  residual  norm.  Because  of  the  compactness  of 
the  integration  in  the  complex  plane,  the  elements  of  the  matrix  can  be  given  in  full  for  a multi- 
component  airfoil  case  and  the  coding  of  a computer  program  becomes  a relatively  simple  matter. 

The  treatment  of  flow  past  airfoils  under  the  constraint  of  wind  tunnel  walls  is  based  on  the 
concept  of  Green’s  function  in  the  complex  plane.  To  a large  extent,  it  is  an  incorporation  of  earlier 
research  results.  References  2-4,  in  the  vortex  panel  method.  The  material  is  organized  in  such  a way 
that  the  method  for  the  calculation  of  flow  past  an  airfoil  in  the  wind  tunnel  is  presented  as  a modifi- 
cation of  the  method  developed  for  an  airfoil  in  free  air,  and  is  easy  to  use  to  assess  the  effects  of  wind 
tunnel  wall  interference.  Although  in  the  context  of  the  present  method  the  introduction  of  the  con- 
cept of  the  angle  of  attack  is  not  necessary  at  all,  it  was  found  to  be  useful  in  the  interpretation  of  the 
wind  tunnel  computations  and  to  show  the  limitations  of  conventional  wall  correction  techniques 
when  applied  to  high  lift  systems. 

2.0  DISTRIBUTION  OF  SOURCES  AND  VORTICES 

Let  D"  be  an  infinite,  connected  region  in  the  complex  plane  bounded  by  several  non- 
intersecting airfoil  contours  Lt , L2 , ...  Ln , oriented  counterclockwise  as  shown  in  Figure  1.  By  L we 
denote  the  union 

n 

L = u Lj 
8 = 1 


, and  by  D+  the  complement  of  D UL.  The  region  D+  thus  consists  of  n simple  closed  interior  regions 
,DJ , D* , ....  D+  bounded  by  the  contours  L, , L2 , ...  Ln . 
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Consider  the  Cauchy  type  integral 


F(z)  = r7 ; dr 

27T1  L S - Z 


where  f is  a point  of  the  contour  L and  f is  a complex  function,  defined  on  L and  continuous  on  each 
of  the  separate  contours  Le . It  is  readily  observed  that  the  Cauchy  type  integral  describes  a sectionally 
analytic  function  F,  which  consists  of  analytic  functions  F“  and  F*  in  the  regions  D-  and  D*,  C = 1,  2, 

. . n.  (The  proof  of  analyticity  consists  in  establishing  the  differentiability  with  respect  to  the  param- 
eter z.)  The  density  function  f determines  the  jump  of  the  function  F across  L,  that  is  the  difference 
of  the  limiting  values  of  functions  F*  and  F~  as  z approaches  the  point  ?0  on  the  contour  Ls : 


Fg(?0 ) ' F'(?o>  " f(fo) 


The  proof  of  Equation  (2.2)  for  a continuous  density  function  f is  given  by  Muskhelishvili,  Reference  5, 
§ 18.  Assuming  that  none  of  the  contours  L8  passes  through  the  point  at  infinity,  the  Cauchy  type 
integral  vanishes  at  infinity: 


F-  (°°)  = 0 (2.2 

The  Cauchy  type  integrals  are  known  to  be  closely  related  to  the  (logarithmic)  potentials 
of  simple  and  double  layers,  see  Reference  5,  § 11.  However,  for  the  purpose  of  modelling  the  flow 
fields  it  appears  more  convenient  to  identify  the  function  F with  the  complex  disturbance  velocity 


F(z)  = U(z)  - iV(z) 


where  U and  V are  the  x and  y-components  of  the  disturbance  velocity.  Since  in  two-dimensional 
potential  flow  the  velocity  field  is  continuous  even  in  the  lifting  case,  there  is  no  need  to  introduce 
discontinuity  lines  (barriers)  in  the  flow  region  D“ . 

It  is  natural  to  ask  now  whether  two  different  density  functions  can  generate  one  and  the 
same  function  F~  in  the  exterior  region  D“ . As  shown  below,  the  answer  is  affirmative. 

Denoting  these  density  functions  fj  and  f2,  Equation  (2.1)  gives  the  corresponding  condition 

— / df  = 0,  zeD" 

2iri  l ? ~ z 


By  the  interior  Cauchy  integral  formula,  Reference  5,  § 24,  the  difference 

*(f)  « f,(f)  ~ f2(f) 

is  then  the  boundary  value  of  the  function 

*(z)  * r^7  / df»  zeD* 

2iri  L i ~ * 


y 
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an?  C°n,tin^°US  in  D+  UL'  In  Practlcal  terms  this  means  that  on  Lfi  the  density 
V ) f *he  ^auchy  tyPe  lntegral  is  determined  by  the  function  F~  to  within  the  boundary  value  ^ (t/ 
of  any  function  (z ) analytic  in  the  interior  region  D* . ^e(i) 

„nn  . . . Con®1<tenng  the  above  degree  of  arbitrariness,  one  can  impose  on  the  density  function  various 

see  Reference  5 1 67  andTsT  rTf^"  °f  ^ b°Utldary  value  Problem  in  Question, 

a,  thTm  " 5’  § 6\  d § 68,  °r  Reference  6,  § 34.  For  the  solutions  to  potential  flow  problems 

bounH^  aPPr°Priate  constraint  appears  the  requirement  that  the  density  function  f(?)  be  the 
thp  r ar^fv  ue  of  a function  f(z),  analytic  in  D and  continuous  in  D“  UL.  Under  such  circumstances 
example^ Reference^^f^l Equation  (2.1),  reduces  to  the  Cauchy  integral,  whose  value  iSS*”™ 


1 

2iri  { 


I f(°°),  zeD 
| -f(z)  + f(°°), 


ream, nap  thJr^u0  fXaminf  the,physical  aspects  of  solutions  based  on  Equation  (2.5),  we  first 
the  Si'  ,ntr0dl,Cine  the  "*le  " *he  "°™al  to  L and 


df  = ie^tn  Idfl  (2.6] 

where  |df  | is  the  contour  length  element.  The  density  function  can  be  expressed  as 

f(f)  = - (a(f)  + (2  7) 

real  functions  o and  y are  called  the  source  and  vortex  densities  respectively.  The  terminol- 
integral  0168  Whe°  Equatlons  <2-6)  and  (2.7)  are  substituted  in  Equation  (2.1).  The  resultant 


PW  * l [ 


contains  singular  terms 


2jt(z  - f)  2jt(z 


%] 


1 __i 

2?r(z  - f ) 2ir(z  - f ) 


thttniTsJenT^  “ d!sturbance  velocities  at  Z due  to  the  unit  strength  source  and 

‘°cated at  Accordingly,  the  integral  of  Equation  (2.8)  will  be  called  the 
distribution  of  sources  and  vortices.  From  the  mathematical  point  of  view,  Equation  (2.8)  is  equiva- 

of  sour ces^^vnrt ’ 1 ' S°i,that  thC  the<>ry  °f  thb  Cauchy  type  integral  can  be  applied  to  the  distribution 
ol  sources  and  vortices  whenever  necessary. 

in  A'  back  Equation  (2.4),  we  can  express  the  component  of  the  disturbance  velocity 

in  the  direction  of  the  unit  vector  exDiiX  las  *y 


U(z)  cos  X + V(z)  sin  X = RelFfzJe**] 


-4- 

Taking  first  X = v and  second  X = v + tt/2,  we  obtain  from  Equations  (2.2)  and  (2.7) 


Re{[Fj(f)  - F-(f)]eiv(f))  = -o(f) 

(2.9) 

Re  { 1 Fg(f ) - F-(f)]eM*«>  + "/2>|  = 7(f) 

(2.10) 

1 

Using  the  above  interpretation,  Equations  (2.9)  and  (2.10)  indicate  that  the  normal  and  tangent  velocity 
jumps  across  L at  f are  given  by  the  values  -o(f)  and  7(f)  respectively. 

Assuming  that  the  free  stream  velocity  vector  is  of  unit  magnitude  and  at  angle  a to  the  real 
axis,  the  normal  and  tangent  components  of  the  resultant  exterior  velocity  can  be  expressed  as  j 

\ 

v(f)  = Re(  [e-i“  + F-(?)]eil'<f>| 

(2.11)  1 

u(f)  = Re  | [e_i“  + F-(f)]eil-'(n^/2][ 

= -Im  |[e-i“  + F-(f)]  e^f  > ) 

(2.12)  | 

The  sign  convention  is  such  that  the  normal  velocity  v is  positive  in  the  direction  of  the  outward  normal 
(with  respect  to  the  airfoil)  and  the  tangent  velocity  u is  positive  in  the  counterclockwise  direction  on 
L. 

Using  Equations  (2.9)  and  (2.10)  we  can  eliminate  F_  from  Equations  (2.11)  and  (2.12): 


v(f)  = Re{[e-i«+  Fj(f)]ei-(f)|+  a(f) 

(2.13) 

u(f)  =-Im|[e-‘“+  Fj(f)]  e^f >|  - 7(f) 

(2.14) 

In  particular,  if 

e-‘“  + Fj(f)  = 0 

(2.15) 

we  obtain 

v(f)  = a(f) 

(2.16) 

u(f)  = -7(f) 

(2.17) 

Equations  (2.16)  and  (2.17)  will  be  satisfied  along  Le  if  Equation  (2.15)  applies  everywhere  on  Le. 

By  the  identity  theorem  then 

e-i°  + Fj(z)  * 0,  zeDj  (2.18) 

/ 


/ 


indicating  that  there  is  no  resultant  flow  in  the  interior  region  Le . 


Conversely,  it  can  be  shown  that  if  either  Equation  (2.16)  or  Equation  (2.17)  is  valid  every- 
where on  Lk  , then  the  other  one  is  valid  too  and  Equation  (2.18)  also  holds  true.  Since  the  establishment 
of  this  fact  plays  a key  role  in  the  development  of  the  present  vortex  panel  method,  the  outline  of  the 
proof  is  presented. 

Assuming  that  Equation  (2.16)  applies  everywhere  on  Lc,  Equation  (2.13)  gives 


e-‘“  + F+(f) 
e-'^f ) 


= 0,  feL„ 


(2.19) 


In  view  of  the  fact  that  both  the  constant  exp  (-  io)  and  the  function  F*  are  analytic  in  D*  Equation 
(2.19)  specifies  a homogeneous  Riemann-Hilbert  problem  for  the  function  <)>*  analytic  in  D+  and 
having  the  boundary  value  6 

*J(f)  * e-io  + FJ(ft,  feLs 

Since  in  going  round  the  contour  Lc  once  in  the  counterclockwise  direction  the  angle  v acquires  the 
increment  2tt,  the  index  of  the  denominator  function,  see  Reference  6,  § 12,  is 

Ind  [e'il,<f)]  = — arg(e~i2,r)  = -1 
2t r 

As  proved  by  Gakhov,  Reference  6,  § 29,  for  a negative  index  the  interior  homogeneous  problem  does 
not  have  a solution.  The  only  possibility  to  satisfy  Equation  (2.19)  is  to  require  that  = 0,  J eLe, 
i.e.  that  Equation  (2.15)  apply  all  over  La.  The  substitution  in  Equation  (2.14)  yields  Equation  (2  17) 
which  completes  the  proof. 

Furthermore,  it  can  be  shown  that  the  distribution  of  sources  and  vortices,  whose  density 
functions  are  subject  to  the  constraints  of  Equations  (2.16)  and  (2.17)  along  the  entire  contour  L,  is 
capable  of  representing  an  arbitrary  function  F~  analytic  in  the  exterior  region  D' . To  prove  this  we 
shall  utilize  the  fact  that  in  this  special  case  the  distribution  of  sources  and  vortices  becomes  the  Cauchy 
integral  in  the  true  sense.  From  Equations  (2.1)  and  (2.18)  we  namely  have 


df  = FUz)  = -e- 


and  consequently,  by  the  exterior  Cauchy  integral  formula,  Equation  (2.5) 


1 f(f) 

_ / jSL  df  , -e-L-w, 


where  f(z)  is  the  analytic  continuation  of  f(f ) into  the  region  D_ , taking  at  the  infinitely  distant  point 


Since  the  exterior  Cauchy  formula  applies  to  an  arbitrary  function  f analytic  in  D~  and  continuous  in 
D"  UL,  then  also 


- 6 - 


is  an  arbitrary  function  analytic  in  D~.  The  density  function  f(f)  is  uniquely  determined  from  the 
limiting  value  of  F~  using  the  relationship 


f(f)  = -e-iQ  - F'(f) 

Having  completed  this  proof,  it  may  be  stated  that  there  is  no  potential  flow  (past  airfoils 
in  free  air)  whose  complex  disturbance  velocity  could  not  be  represented  by  the  distribution  of  sources 
and  vortices,  satisfying  Equations  (2.16)  and  (2.17). 

Besides  expressing  the  jumps  in  normal  and  tangent  velocities  according  to  Equations  (2.9) 
and  (2.10),  the  source  and  vortex  densities  also  acquire  a distinct  physical  meaning  when  integrated 
along  L,  namely  that  of  the  discharge  and  of  the  circulation. 

Since  the  functions  F+  are  analytic  in  the  simple  connected  interior  regions  D^,  they  are 
subject  to  the  Cauchy  integral  theorem 


/ Fj(f ) df  = 0,  8 = 1 n 

Ls 

Thus,  by  the  integration  of  the  jump  condition.  Equation  (2.2),  and  with  the  help  of  Equations  (2.6) 
and  (2.7)  we  have 


/ F-«-)dr 

L 


/ f(f ) df 

L 


/ 7(f)ldf|  + i / a(f)|df| 
L L 


(2.20) 


From  Equations  (2.11),  (2.12),  and  (2.6)  also 


/ F-tf)dr  = /u(f)|df|  + i / v(f)|df| 

L L L 


= -T  + iQ 


(2.21) 


where 


Q = / v(f)ldf I = / o(?)ld?l 

L L 


(2.22) 


is  the  discharge  (efflux)  through  L and 


r = - / u(r)|df|  = / 7(f)ldfl 

L L 


(2.23) 


is  the  circulation  around  L.  Of  course,  in  the  case  that  Equations  (2.16)  and  (2.17)  apply,  Equations 
(2.22)  and  (2.23)  are  immediately  apparent. 


1, 


The  discharge  and  the  circulation  play  an  important  role  in  the  far  field  approximation  of 
the  complex  disturbance  velocity . Assuming  that  the  airfoil  is  located  near  the  point  z() , we  can  expand 
F"  for  large  |z  - zu  | in  the  Laurent  series 


where 


F (z)  = £ cm  (z  - z0)' 

m = 1 


= — / «■-  z0r_l  F-(f)  dr 


The  constant  term  cD  is  missing  because  of  Equation  (2.3).  Using  Equation  (2.21)  we  have 


F (Z)  = q— + r + 

2tt(z  - zG ) 2tt(z  - z0 ) 


(z  - ZQ  ) ->  oo 


(2.24) 


Hence,  far  from  the  airfoil  the  complex  disturbance  velocity  can  be  approximated  by  a point  source 
of  strength  Q and  by  a point  vortex  of  strength  f,  located  at  zc . 

3.0  AIRFOIL  BOUNDARY  VALUE  PROBLEM 

Assuming  the  free  stream  velocity  to  be  of  unit  magnitude  and  at  angle  a to  the  real  axis, 
the  airfoil  boundary  condition  can  be  written,  in  accordance  with  Equation  (2.11) 

Re  | [ e~ ia  + F-tf))^?)}  = v(r)  (3.1) 

The  normal  component  v of  velocity  is  regarded  as  a given  function  of  L.  If  the  flow  is  tangent  to  the 
contour,  then  v(r)  = 0.  The  use  of  nonzero  values  of  v serves  the  purpose  of  approximating  the  dis- 
placement effect  of  thin  boundary  layers  in  terms  of  flow  through  the  airfoil  contour.  Based  on  the 
concept  of  the  displacement  flux  by  Preston,  Reference  7,  the  normal  component  of  velocity  is  set 
equal  to  the  stream  wise  derivative  of  the  product  of  the  local  potential  flow  velocity  and  the  displace- 
ment thickness  of  the  boundary  layer;  for  more  details  see  Reference  8.  Compared  to  the  classical 
approach  by  Prandtl,  where  the  displacement  thickness  is  added  to  the  airfoil  thickness  to  produce  a 
modified  shape,  the  Preston  concept  has  the  advantage  of  preserving  airfoil  geometry  in  the  course  of 
potential  flow-boundary  layer  iterations. 

The  condition  of  finite  trailing  edge  velocity  — the  Kutta-Joukowki  condition  — can  be 
derived  in  terms  of  the  source  and  vortex  densities  using  the  properties  of  the  Cauchy  type  integral, 
Equation  (2.1).  From  the  jump  condition  of  Equation  (2.2)  it  is  apparent  that  for  the  establishment 
of  the  finite  velocity  field  the  density  f has  to  be  finite  everywhere  on  L.  Since  a discontinuity  in  f 
gives  rise  to  a logarithmic  singularity  of  the  Cauchy  type  integral,  see  Reference  6,  § 8,  we  require 
that  the  density  f be  continuous  on  each  of  the  airfoil  contours  Le , including  the  trailing  edge  point. 
Denoting  by  subscripts  U and  L the  upper  and  lower  trailing  edge  values,  we  thus  require  that 


fU  = fL 


Applying  Equation  (2.7),  we  have 


(®u  + iyvle'^u  = (oL  + i7L)  e-,*'L 
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or 


°L  + iTl  = ~(Ov  + i7u)eie 

(3-3) 

where 

0 = vL  - vv  - 1 T 

(3.4) 

is  the  trailing  edge  angle,  see  Figure  2.  Thus,  Equation  (3.3)  may  be  regarded  as  the  Kutta-Joukowski 
condition  expressed  in  terms  of  the  source  and  vortex  densities. 

Independent  conditions  for  the  source  and  vortex  densities  are  obtained  from  Equation  (3.3) 
in  the  following  special  cases: 

1.  blunt  trailing  edge,  6 = it 

°l  = °u 

(3.5) 

TT  = T'u 

(3.6) 

2.  cusped  trailing  edge,  0 = 0 

' °L  = ~°v 

(3.7) 

Tl  = -7u 

(3.8) 

However,  if  0 < 0 < n,  Equation  (3.3)  gives 

°L  + °u  cos  ® 

7v  ~ ■ a 

sin  0 

(3.9) 

°u  + °L  cos  ® 

7L  " . a 

sin  6 

(3.10) 

indicating  that  the  trailing  edge  values  of  the  source  and  vortex  densities  are  interdependent. 


Geometrically,  Equations  (3.9)  and  (3.10)  express  the  condition  of  equal  velocities  at  the 
upper  and  lower  trailing  edge  points.  This  can  easily  be  established  by  considering  the  velocity  triangles 
in  Figure  2.  We  have 


"uu  vu  VL  UL 
cos  flu  sin  0^  sin  0L  cos  0L 


where 
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Eliminating  and  0L  and  expressing  u,j  and  uL  in  terms  of  Vy , vL , and  0 we  arrive  at  the  equations 


uu 


VL  +-»u  COS  0 
sin  6 


(3.11) 


UL 


VU  + VL  cos  ® 
sin  0 


(3.12) 


Substituting  for  the  normal  and  tangent  velocities  from  Equations  (2.13)  and  (2.14)  and  using  Equation 
(3.4),  we  again  obtain  Equations  (3.9)  and  (3.10).  In  the  case  that  Equations  (2.16)  and  (2.17)  apply, 
the  operation  is  particularly  simple;  this  way  of  deriving  the  trailing  edge  conditions  was  first  used  by 
Dvorak  and  Woodward,  Reference  9.  Hence,  with  respect  to  the  trailing  edge  velocity  the  concepts  of 
finiteness  and  single-valuedness  are  found  to  be  equivalent. 

Going  back  to  Equations  (3.11)  and  (3.12),  another  important  observation  can  be  made.  If 
the  normal  components  of  velocity  are  zero,  then  also  the  tangent  components  vanish  and  the  trailing 
edge  is  inevitably  a stagnation  point.  The  only  exception  to  this  rule  is  a cusped  trailing  edge,  where 
Equations  (3.11)  and  (3.12)  do  not  apply.  For  a cusped  trailing  edge,  6 = 0,  we  obtain  from  Equations 
(3.7)  and  (3.8)  the  velocity  relationships 


vL  = "vu 
Ul  = "Uu 


which  do  not  imply  the  existence  of  a stagnation  point  if  vL  * Vu  = 0.  The  best  known  example  is  the 
Joukowski  airfoil.  For  a blunt  trailing  edge,  6 = jt,  we  obtain  from  Equations  (3.5)  and  (3.6)  the 
relationships 


VL  = VU 


uL  = uu  = 0 (by  definition) 
so  that  again  a stagnation  point  is  obtained  if  vL  = vy  = 0. 

In  conclusion  of  this  section  a few  general  comments  will  be  made  concerning  the  methods 
of  solution  of  the  airfoil  boundary  value  problem  specified  by  Equations  (3.1)  and  (2.3).  By  substitu- 
ting in  Equation  (3.1)  from  Equation  (2.8)  we  obtain  a singular  integral  equation  which  contains  two 
unknown  real  functions  a and  7.  From  the  discussion  in  Section  2.0  it  became  already  apparent  that 
the  most  natural  approach  is  to  eliminate  the  source  density  o by  setting  it  equal  to  the  prescribed 
normal  component  of  velocity.  Equation  (2.16).  It  was  shown  that  this  type  of  constraint  is  allowable 
since  the  integral  of  Equation  (2.8)  is  still  capable  of  representing  an  arbitrary  function  F~  analytic  in 
the  flow  region  D“ . In  view  of  the  fact  that  the  vortex  density  7 is  the  function  to  be  evaluated,  we 
speak  of  the  vortex  distribution  method.  Its  important  feature  is  that  the  Kutta-Joukowski  condition 
for  a finite-angle  trailing  edge  is  satisfied  by  simply  fixing  the  trailing  edge  values  of  7 from  the  known 
values  of  a according  to  Equations  (3.9)  and  (3.10).  As  shown  in  Section  2.0,  the  evaluated  distribution 
of  7 gives  directly  the  distribution  of  the  tangent  component  of  velocity.  Equation  (2.17). 

In  contrast  to  the  above  approach,  the  source  distribution  method  utilizes  the  vortex  density 
as  a prescribed  function,  and  the  source  density  is  the  quantity  to  be  solved  for.  Apart  from  the  re- 
quirements of  continuity  on  smoth  portions  of  the  contour  and  of  a sufficient  degree  of  freedom  to 
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satisfy  the  trailing  edge  condition,  the  vortex  density  can  be  prescribed  arbitrarily*.  Since  neither 
of  the  trailing  edge  values  of  o and  7 are  known  a priori,  the  Kutta-Joukowski  condition  for  a finite- 
angle  trailing  edge  is  more  difficult  to  satisfy.  The  most  frequently  used  technique,  also  known  as  the 
Douglas-Neu'mann  method,  References  1 and  10,  satisfies  the  trailing  edge  condition  in  terms  of  the 
linear  superposition  of  the  “basic  flows”:,  the  nonlifting  flows  due  to  uniform  streams  at  0 and  90 
degree  incidences,  the  circulatory  flows  due  to  a prescribed  distribution  of  7 on  each  airfoil  component, 
and  the  displacement  flow  due  to  the  prescribed  normal  component  of  velocity  along  the  airfoil  con- 
tour. Using  this  approach,  Equations  (2.16)  and  2.17)  are  not  satisfied  and  therefore  the  tangent 
component  of  velocity  has  to  be  obtained  by  evaluating  the  integral  of  Equation  (2.8).  A complex 
plane  version  of  this  method,  using  flat  panels  with  constant  source  and  vortex  densities,  is  described 
in  Reference  1 and  no  attempt  has  been  made  to  duplicate  this  effort  here.  The  vortex  distribution 
method  is  not  only  more  elegant  and  straightforward  but,  as  the  comparisons  of  computer  results  in 
Section  7.0  reveal,  also  more  accurate. 


4.0  VORTEX  PANEL  METHOD 

Each  of  the  n airfoil  components  is  approximated  by  an  inscribed  polygon  as  shown  on  the 
C-th  component  in  Figure  3.  The  corner  points  are  numbered  counterclockwise,  starting  with  the 
(upper)  trailing  edge  point,  j = m(8  - 1)  + 1,  and  ending  again  with  the  (lower)  trailing  edge  point, 
j = m(8).  Since  on  the  first  airfoil  component  the  number  of  the  upper  trailing  edge  is  1,  we  define 
m(0)  = 0.  The  last  point  is  the  lower  trailing  edge  point  of  the  n-th  component;  m(n)  is  therefore  the 
total  number  of  corner  points.  The  density  functions  a and  7 are  assumed  to  vary  linearly  over  each 
panel  (line  segment),  be  continuous  across  the  comer  points  on  the  smooth  portions  of  the  airfoil 
contour,  j = m(8  - 1)  + 2,  m(8  - 1)  + 3, . . .,  m(8)  - 1,  and  obey  the  trailing  edge  conditions,  Equations 
(3.9)  and  (3.10),  for  points  j = m(8  - 1)  + 1 and  m(8).  The  trailing  edge  angle  is  supposed  to  satisfy 
0 < 0S  < IT. 

Assuming  that  all  comer  values  Oj  are  known  from  the  prescribed  distribution  of  the  normal 
velocity,  this  leaves  m(8)  - m(8  - 1)  - 2 unknown  values  y-.  of  the  vortex  density  on  the  8-th  airfoil 
component  and  m(n)  - 2n  unknown  values  7j  for  the  whole  airfoil.  They  can  be  found  by  satisfying 
the  boundary  condition.  Equation  (3.1),  at  panel  midpoints.  Since  there  are  m(n)  - n panels,  the 
resultant  system  of  linear  algebraic  equations  is  overdetermined  by  n equations.  By  solving  the  over- 
determined system,  we  have  the  boundary  conditions  satisfied  for  all  panels,  but  with  some  error 
which  is  subject  to  minimization.  This  idea,  suggested  to  the  author  by  his  colleague  D.J.  Jones, 
appears  to  have  distinct  advantages  over  sacrificing  n selected  panels  by  dropping  their  boundary 
conditions  and  concentrating  the  error  there. 

Using  Equation  (2.1)  we  can  express  the  complex  disturbance  velocity  at  the  flowfield 

point  z as 


n m(C)-l 

F-(z)  = 2:  r A;F-(Z)  (4.1) 

e = 1 j = m(S-l )+  1 


The  Kutta-Joukowski  condition  indirectly  implies  that  the  distribution  of  the  vortex  density  has 
to  provide  the  proper  amount  of  circulation  around  each  airfoil  component.  Apart  from  this, 
there  are  no  additional  conditions  to  consider  when  prescribing  the  vortex  density.  As  in 
Reference  6,  § 34,  the  proof  can  be  provided  from  the  conditions  of  the  solubility  of  the 
corresponding  Riemann-Hilbert  problems  in  the  interior  regions  Dj. 


where 


Aj  F~  (z)  = 


2iri 


rj+ 1 
/ 


f(f) 

?-  z 


df 


(4.2) 


is  the  contribution  of  the  j-th  panel.  Taking  into  account  the  linear  dependence 


m = fj 


fj+i 


fj+ 1 - 


(f-fj) 


we  obtain 


fj+i-fj 


(?-  Z)  + fj+1 


z~fj 
fj+l  ” fj 


z fj+l 

fj+1  ^ 


J_  fi+1  J(f)_  = Vi  ~ fj  + JL  l\  z~  fj 

2*i  f-z  2rri  + 2rri  [_  1+1  fj+1  - fj 


fj 


z " fj+l 
fj+ 1 “ fj  . 


log 


fj+l  z 
fj  - z 


(4.3) 


From  Equation  (2.7)  we  have  for  the  j-th  panel 


fj  = ~(Oj  + »7j)  e"il'j 


(4.4) 


fj+l  = ~ (°j+ 1 + i7j+i ) e_il'j 


(4.5) 


where,  according  to  Equation  (2.6) 


e^j 


. fj+l  ~ fj 
'lfj+l  -fjl 


(4.6) 


Since  AjF"  is  supposed  to  be  discontinuous  as  z crosses  the  segment  f j » f j +- 1 and  vanish  at 
z = °°,  we  select  the  negative  semiaxis  as  the  branch  line  of  the  logarithmic  function  and  require  that* 


log  (1  + Oi)  = 0 


As  a result 


fj+l  1 

fj“z 


+ ir 


As  a rule,  computer  subprograms  for  the  complex  logarithmic  function  are  written  to  evaluate  the 
“principal  value”,  satisfying  these  requirements. 
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where  - n < r < it  is  the  angle  obtained  by  the  rotation  of  the  vector  f ) - z into  the  direction  of  the 
vector  jj+1  - z,  see  Figure  3.  In  the  limit  as  the  point  z approaches  any  point  of  the  j-th  panel  from 
the  flowfield  side,  r tends  to  the  value  -7r.  Thus  the  complex  disturbance  velocity  induced  by  the 
panel  at  its  own  midpoint  is 


A TT—  ^ +^)+'\  ^ 1 

i,r\— 2 )"  ~2ri~  ' 4 <'J*'  +,i> 


Denoting  the  midpoint  of  the  k-th  panel 


zk  g tfk  + fk  + i ) 


and  using  Equations  (4.1)  to  (4.7),  we  obtain 


where 


m(n) 

F (zk)  = 2 Ckj  (0j+i7j) 

j = l 


Ck,j  " Kk,j  + Lk,j 


1 zk  ?j  + l ?j+l  ~ zk 

— e—j  1 + log  , j * k,  m(S) 

2m  L ?i+i  fj " zk  J 


. j = k 


(4.10) 


(4.11) 


, j = m(2) 


1 , T fj~zkl 

|(i ■i'k*1 


(4.12) 


' * j = m(8-  1)+  1 

Using  Equation  (2.16)  and  applying  the  airfoil  boundary  condition,  Equation  (3.1),  at  panel 
midpoints  zk , we  obtain 


Re  J Ie-‘“  + F^(zk)]  ei*’k  j - j (ok  + <7k+1) 


(4.13) 


or,  with  the  help  of  Equation  (4.9) 


»>(n)  r m (n)  "II 

2 Im  (e'^kCk  j)7j  = Re  ^eil'k(e"i * *“  + 2 Ck  j 0j)J  -- (ok  + ak  + , ),  k* 

j = 1 j = 1 2 

Introducing  the  index  K = k + £,  this  equation  can  be  written  in  the  matrix  form 


AKj7j  = bK , K = 1,  2, . . m(n)  + n 


j = 1,  2, . . m(n) 


where 


(4.14) 


Ak+B j = Im  (e'^k  Ckj),  ktm({-  l),m(£) 


(4.15) 


[m  (n ) -i  X 

e"'k(e-1“  + 2 Ck>jOj)  --(ok+ak  + 1),  k 4 m(C  - 1),  m(C)  (4.16) 

j = l J 4 


The  elements  of  rows  m(C- 1)  + £ and  m(£)  + £ are  obtained  from  the  trailing  edge  conditions,  Equations 
(3.9)  and  (3.10): 


AmCs-n+aj  = 


1,  j = m(£-l)  + 1 


0,  otherwise 


(4.17) 


°m  (8)  + °m  (8-  1 )+  1 C°S08 


(8-1  )+8 


(4.18) 


Vn(8)+e,j 


1,  j - m(fi) 


0,  otherwise 


(4.19) 


3m(S)+e,j 


7m(8-l)+l  + am(8)c°8^e 
sin  0B 


(4.20) 


Utilizing  the  rotation  of  unit  trailing  edge  vectors,  the  cosine  and  the  sine  of  the  trailing  edge  angle 
6 1 are  obtained  as  the  real  and  imaginary  parts  of  the  expression 


cos  On  + i sin  0B 


?m(8)-l  " ?m(8)  Ifm  (8-  I )+2  " ?m(8-l)+l  I 
lfm(8)-l  -fm(8)l  fm  (8- 1 )+2  " (8- 1 )+l 


(4.21) 
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To  obtain  a least  squares  solution  to  the  overdetermined  system  described  by  Equation  (4.14), 
we  first  premultiply  both  sides  by  the  transpose  matrix  Aj  K : 


<Aj,K  AKJ)7j  = Aj  K bK 


(4.22) 


Assuming  that  the  panels  do  not  overlap  or  cross,  that  is  have  no  common  points  other  than  the  corners, 
the  [m(n)  + n]  X m(n)  matrix  AK  j is  of  rank  m(n).  Then  the  square  m(n)  X m(n)  matrix ( Aj  K AK  j ) is 
nonsingular  and  of  rank  m(n).  Consequently,  the  system  of  Equation  (4.22)  possesses  a unique  solution 
which  may  be  found  by  the  LU  factorization  or  other  efficient  methods  for  solving  systems  of  linear 
algebraic  equations.  Using  the  fact  that  the  matrix  (Aj  K AK  j ) is  symmetric  positive-definite,  one  can 
use  with  great  advantage  Crout’s  factorization  scheme,  Reference  11,  § 2.3.4. 1,  selecting  the  pivots  in 
order  along  the  diagonal,  as  suggested  in  Reference  12. 

Denoting  by 

rK  = bK  " ak  ,j  Tj 

the  residual  of  the  original  system,  it  can  be  shown  see  for  example  Reference  13,  § 5.6,  that  the 
solution  of  Equation  (4.22)  minimizes  the  sum 


m (n)+n 


Of  course,  there  are  other  methods  for  solving  the  overdetermined  system  of  Equation  (4.14), 
such  as  the  minimization  of  the  Lj  norm  of  the  residuals. 


m (n)+n 

2 |rK  I 

K = 1 

or  the  minimization  of  the  maximum  residual  (L„  or  Chebyshev  norm) 


max  |rK  | 
K 


These  solutions  can  also  be  very  efficient  when  based  on  linear  programming  techniques,  see  References 
14  to  17.  As  demonstrated  on  examples  in  Section  7.0,  from  the  aerodynamic  point  of  view  they  give 
almost  identical  results  as  the  least  squares  solution.  However,  they  can  be  of  great  value  especially  if 
the  system  (4.14)  is  ill-conditioned  due  to  the  occurence  of  nearly  coinciding  panels,  since  the  condi- 
tion of  Equation  (4.22),  on  which  the  least  squares  solution  is  based,  is  worse  than  that  of  Equation 
(4.14),  see  Reference  13,  § 8.8.  Otherwise,  the  least  squares  solution  offers  better  economy  when 
repeated  computations  for  the  same  airfoil  configuration  are  required,  such  as  in  iterative  coupling  of 
the  potential  flow  and  boundary  layer  computations.  The  advantage  of  the  least  squares  solution  is 
that  it  is  possible  to  generate  the  inverse  of  the  matrix  AjK  AK  j so  that  the  subsequent  computations 
are  reduced  to  the  matrix  multiplication  of  the  inverse  by  the  new  right  hand  side  Aj  K bK . 

Once  the  values  7j  of  the  vortex  density  have  been  obtained,  the  comer  values  of  tangent 
velocity  are  known  from  Equation  (2.17).  Using  the  Bernoulli  theorem,  the  pressure  coefficients  are 


evaluated  as 


Cp(fj)  = 1 - (o2  + y2)  (4.23) 

Admittedly,  Equations  (4.1)  to  (4.6)  indicate  at  the  same  time  that  the  complex  disturbance 
velocity  is  singular  at  the  corner  points  fj  of  the  polygon,  with  the  exception  of  the  trailing  edge  where 
the  Kutta-Joukowski  condition  is  satisfied.  However,  the  trailing  edge  is  real  whereas  other  comers  lie 
on  smooth  portions  of  the  original  airfoil  contour.  If  we  accept  the  panel  method  just  as  a numerical 
technique  which  allows  us  to  discretize  the  continuous  mathematical  model  in  order  to  evaluate  the 
flow  quantities  at  the  points  of  a smooth  contour  — and  not  at  corner  points  of  an  approximation 
polygon  — the  evaluation  of  the  pressure  coefficient  from  Equation  (4.23)  is  justifiable.  Comparisons 
with  exact  flow  solutions  fully  confirm  this  philosophy.  The  advantage  of  the  present  approach  is  that 
the  pressure  coefficients  are  evaluated  directly  at  the  contour  input  points  whereas  most  panel  methods 
give  output  at  panel  midpoints,  which  in  fact  are  not  points  of  the  contour. 


5.0  DISTRIBUTION  OF  SOURCES  AND  VORTICES  MODIFIED  TO  INCLUDE 
THE  WALL  INTERFERENCE  EFFECT 

It  is  assumed  that  the  flow  region  is  bounded  from  the  outside  by  two  parallel,  infinite  wind 
tunnel  walls  at  y = h/2  and  -h/2,  see  Figure  4,  and  that  the  walls  behave  in  such  a way  that  the  normal 
component  of  velocity  is  proportional  to  the  pressure  drop  across  the  wall  (Darcy  law).  We  have  mainly 
the  perforated  walls  in  mind,  but  indications  are  that  even  slotted  walls  behave  to  some  extent  like  this, 
Reference  18. 

If  the  pressure  outside  the  walls  (the  plenum  pressure)  is  equal  to  the  static  pressure  of  the 
wind  tunnel  stream  infinitely  far  upstream,  then  the  normal  component  of  velocity  is  proportional 
directly  to  the  wall  pressure  coefficient.  In  the  wind  tunnel,  the  flow  infinitely  far  upstream  is  sup- 
posed to  be  at  angle  a = 0 with  respect  to  the  x-axis,  so  that  in  terms  of  complex  disturbance  velocity 
the  normal  component  of  velocity  is  -Im  F~  on  the  upper  wall  and  Im  F~  on  the  lower  wall.  Assuming 
further  that  the  airfoil  induces  small  pressure  disturbances  at  the  walls,  the  pressure  coefficient  can  be 
linearized  and  expressed  as  - 2 Re  F" . Summarizing  this,  we  obtain  the  following  boundary  conditions 
for  upper  and  lower  walls 


Py  Re  F"(z)  - 

Im  F-(z)  = 

0, 

z = x + i h/2 

(5.1) 

PL  ReF-(z)  + 

ImF-(z)  = 

0, 

z = x - i h/2 

(5.2) 

The  positive  constants  Py  and  PL , called  the  porosity  parameters  or  resistance  factors  for  upper  and 
lower  walls,  have  to  be  determined  experimentally  for  given  walls  and  specific  flow  conditions;  either 
directly.  Reference  18,  or  indirectly,  Reference  4.  The  solutions  based  on  Equations  (5.1)  and  (5.2) 
contain,  as  limiting  cases,  those  for  a solid  wall  wind  tunnel  (Py  = PL  = 0)  and  for  an  open  jet  wind 
tunnel  (Py  ,PL-*  °°).  The  reason  for  using  different  porosity  parameters  for  upper  and  lower  walls  is 
that  the  walls  may  actually  be  of  different  open  area  ratios  or,  even  if  identical,  they  may  display  a 
different  resistance  to  outflow  than  to  inflow,  see  for  example  Reference  4.  Under  the  latter  circum- 
stance, Equations  (5.1)  and  (5.2)  are  not  suited  for  mixed  inflow-outflow  conditions,  but  apply  well 
if  along  each  wall  either  the  inflow  or  outflow  conditions  prevail,  as  is  the  case  of  high  lift  systems 
(inflow  on  the  suction  side  and  outflow  on  the  pressure  side  of  the  airfoil).  * 

Since  Equations  (5.1)  and  (5.2)  are  linear  and  homogeneous  with  respect  to  the  real  and 
imaginary  parts  of  F- , we  may  use  the  concept  of  a Green’s  function  to  express  the  complex  disturb- 
ance velocity  of  the  wind  tunnel  flow  past  a multi-component  airfoil  as  a line  integral  along  the  airfoil 
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contour.  In  the  context  of  Green’s  theorems  and  eigenfunction  expansions,  the  basic  idea  of  this 
approach  was  first  given  by  Kacprzynski,  Reference  19. 

Assuming  that  the  airfoil  is  positioned  inside  the  wind  tunnel, 


Hmfl  < h/2,  feL 
we  can  write,  in  analogy  to  Equation  (2.8) 

F(z)  = / [a(f)G0(z,f)  + 7(f)G7(z.f)]  Idfl 

L 

where 


G0(z,n  = 

1 

+ H0(Z,n 

to 

N 

G7(z,f)  = 

i 

+ H7(2,n 

2ir(z-  f) 

(5.3) 

(5.4) 

(5.5) 


are  the  Green’s  functions  (influence  functions). 

The  first  right  hand  terms,  generally  called  the  fundamental  solutions,  are  singular  at  z = f 
and  vanish  at  z = °°.  They  are  recognized  as  complex  disturbance  velocities  at  z due  to  a point  source 
and  a point  vortex  at  f , as  if  no  walls  were  present.  The  functions  H0  and  H7  are  analytic  in  the  strip 
|Imz|  < h/2  and  such  that  both  G0  and  Gy  satisfy  the  wall  boundary  conditions  of  Equations  (5.1) 
and  (5.2).  In  terms  of  a two-argument  Green’s  function  G (=  G„  or  G7 ) these  conditions  read 


?u  Re  G(z,f ) - Im  G(z,f)  = 0,  z = x + i h/2 


(5.6) 


1 

1 


i 


I'L  Re  G(z,f)  + Im  G(z,f ) =0,  z = x - i h/2  (5.7) 


for  any  |Imf|  < h/2.  Consequently,  the  functions  Ha  and  Hr  can  be  interpreted  as  complex  disturbance 
velocities  at  z induced  by  the  wind  tunnel  walls  as  a response  to  the  presence  of  a point  source  and  a 
point  vortex  at  f,  respectively.  Using  Equations  (5.4)  to  (5.7)  it  can  be  shown  that  the  integral  of 
Equation  (5.3)  satisfies  the  wall  boundary  conditions.  Equations  (5.1)  and  (5.2),  for  arbitrary  density 
functions  a and  y. 


Solving  the  boundary  value  problem  of  Equations  (5.6)  and  (5.7)  for  the  Green’s  functions 
defined  by  Equations  (5.4)  and  (5.5)  by  the  image  method  of  Ebihara,  Reference  20,  the  analytic  parts 
are  found  to  be.  Reference  3: 


H (z.t)  = R (z.t)  + E (z.t)  + C 15.81 
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where 


B(z,f ) = 


2h 


'U  + ,L  z-f 

IT 

e 2 h 
Z-f 

w 

e h - 1 


. 1 

27r(z-f) 


(5.10) 


E(z,f)  = 


'U*1!-  z-f 

_L  e 2 h 

2h  z-f 

ir 

e h +1 


(5.11) 


C = 


0 , otherwise 


' 

(5.12) 


In  Equation  (5.11),  f denotes  the  complex  conjugate  of  f . The  porosity  parameters  P,,  and  PL  enter 
the  above  formulas  through  the  functions 


" 


*L  = 


which  vary  on  the  interval  < 0,1 ). 

Further  details  concerning  the  derivation  are  not  given  here,  but  the  results  can  easily  be 
verified  by  direct  substitution  in  Equations  (5.6)  and  (5.7).  From  Equations  (5.6)  and  (5.7)  it  is  also 
apparent  that  the  constant  C can  acquire  an  arbitrary  real  value  if  Pu  = Pl  = 0;  otherwise  it  must  be 
zero.  The  value  C = 1/h  selected  for  the  solid  wall  wind  tunnel,  Py  = Pl  = 0,  satisfies  the  additional 
condition 

lim  Gg  (z,f)  = 0 
Re  z -*  -°° 

by  which  the  undisturbed  flow  conditions  are  established  infinitely  far  pstream  (but  not  necessarily 
far  downstream).  If  at  least  one  of  the  porosity  factors  is  non-zero,  then 


— atan  P, 

7T 


2 

— atan  PL 
7T 


(5.13) 


(5.14) 


lim  B(z,f)  = lim  E(z,f)  = 0 
Rez-*±°°  Re  z ->■  ± «= 

so  that  undisturbed  flow  conditions  are  automatically  established  both  upstream  and  downstream. 
This  basic  difference  between  the  solid  wall  and  the  porous  wall  wind  tunnels  is  due  to  the  simplifying 
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assumption  that  the  walls  are  of  infinite  streamwise  extent,  see  Reference  21.  In  the  porous  wall  wind 
tunnel  the  flux  due  to  a source  can  escape  through  the  walls,  no  matter  how  little  porous  they  are,  and 
the  undisturbed  flow  conditions  are  established  infinitely  far  upstream  and  downstream.  For  the  solid 
wall  wind  tunnel  this  is  not  true  and  the  presence  of  a source  is  thus  felt  over  the  entire  wind  tunnel 
length.  The  only  exception  is  if  the  discharge  Q,  Equation  (2.22),  of  the  airfoil  is  zero,  so  that  the 
wall  interference  contributions  of  the  source  terms  cancel  out  far  upstream  and  downstream.  The 
awareness  of  the  above  properties  of  the  wind  tunnel  flow  solutions  is  very  important  for  the  correct 
interpretation  of  the  results. 

It  is  easily  verified  that  the  functions  H0  and  H7  have  no  singularities  other  than  simple  poles 
outside  the  wind  tunnel.  Their  locations 


z = ? ± i2kh  , k = 1,2,  3, ... 


z = ? ± i(2k - l)h,  k = 1,2,3,... 

‘ 

are  known  as  “images”  of  f . At  z = f the  function  B has  a removable  singularity,  so  that  the  functions 
Ha  and  H are  analytic  in  the  entire  wind  tunnel  interior.  In  order  to  calculate  their  values  in  the 
vicinity  of  point  z = f , we  use  the  power  series  expansion,  Reference  2: 


S(z,r)  = — 

2h 


iz-  n 

h 


< 2 


(5.15) 


where  Bm  (t)  denotes  the  Bernoulli  polynomial 


In  particular,  at  z = f we  obtain 


B(z,z) 


2h 


J_  /lu  +tL  1 \ 

2h  \ 2 2/ 


(5.18) 


If  the  wind  tunnel  height  h becomes  infinitely  large,  the  functions  H0  and  H7  vanish,  and 
Equation  (5.3)  reduces  to  Equation  (2.8),  valid  for  the  free  air  case. 

Since  the  functions  H0  and  H7  are  analytic  in  the  entire  wind  tunnel  interior,  they  must  also 
be  continuous  across  the  airfoil  contour  L and  hence  cannot  affect  the  jump  condition  described  by 
Equation  (2.2).  From  this  it  already  follows  that  Equations  (2.9)  to  (2.19)  and  conditions  under  which 
they  were  derived  remain  valid.  Consequently,  in  the  vortex  distribution  method  the  evaluation  of  the 
vortex  density  7 from  the  airfoil  boundary  condition,  Equations  (3.1)  and  (2.16),  and  from  Equations 
(5.3)  to  (5.5)  remains  essentially  the  same  as  for  the  free  air  case.  At  this  stage  attention  is  no  longer 
paid  to  the  wind  tunnel  boundary  conditions,  Equations  (5.1)  and  (5.2),  since  they  are  automatically 
satisfied  by  the  properties  of  the  above  Green’s  functions  for  any  distribution  of  sources  and  vortices. 
As  in  the  free  air  case,  the  tangent  component  of  velocity  is  obtained  from  Equation  (2.17)  and  the 
pressure  coefficient  follows  from  Equation  (4.23). 

If  the  airfoil  chord  to  wind  tunnel  height  ratio  c/h  is  small,  say  1/4  or  less,  the  wall  inter- 
ference effect  in  potential  flow  can  be  approximated  by  a single  angle  of  attack  correction.  Using  the 
analytic  part  of  the  integral  of  Equation  (5.3),  the  angle  of  attack  correction  can  be  expressed  as 


Aa  = -Imf  [o(f)H0(z,f)  + 7(f)H7(z,r)]  |d?|  (5.19) 

where  z is  some  reference  point  inside  the  airfoil,  e.g.  the  quarter  chord  point.  In  order  to  obtain  a 
formula  involving  only  the  integrals  of  a and  7,  we  expand  the  functions  Hc  and  H7  in  Taylor  expan- 
sions about  the  point  z: 


H(z,f)  = H(z,z)  + 


dH(z,n 

d? 


(f  - z)  + . . . 
r** 


Assuming  that  the  angle  of  attack  correction  does  not  vary  significantly  along  the  airfoil,  we  restrict 
ourselves  only  to  the  constant  terms.  Thus 


Aa  = -Im{H0(z,z)/o(r)ldf|  + H7(z,z)  / 7(f)ldf l} 

L L 

= - Q Im  H0(z,z)  - rimH7(z,z) 


where  Q and  T,  Equations  (2.22)  and  (2.23),  are  the  values  of  the  discharge  and  the  circulation,  respec- 
tively. Of  course,  the  above  expression  can  also  be  obtained  from  the  far  field  approximation  of  F~ 
in  free  air,  Equation  (2.24),  employing  the  source  and  the  vortex  strengths,  Q and  f,  as  the  factors  of 
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the  analytic  parts  of  Green’s  functions,  Equations  (5.4)  and  (5.5).  The  values  H0  (z,z)  and  H (z,z)  can 
be  evaluated  from  Equations  (5.8),  (5.9),  (5.11),  and  (5.18).  In  particular,  if  the  airfoil  is  located  near 
the  tunnel  axis,  then  z = 0 and 


Q . / tu  - tlA  P r ' / - h\ 

Aa‘Th‘my—)-  J 


If  the  flow  is  tangent  to  the  airfoil,  then  Q = 0 and  according  to  the  Joukowski 


- l]  (5.20) 

theorem  approximately 


r = 


1 

— c 

2 


CL 


where  CL  is  the  lift  coefficient  of  the  airfoil.  As  a result 


1 c tu  + fcL  1 r / tu  " 4l  \ -i ) 

A a C.  + — cos  ( 7r ) - 1 } 

4 h L 2 2 L V 2 / J| 


(5.21) 


which  in  the  case  of  equal  porosity  parameters,  Py  = PL  = P,  that  is  ty  = tL  = t,  reduces  to  the  familiar 
formula 


A a 


]_ 

4 


iCl6 


L°0 


(5.22) 


The  expression 


— t atan  P 

4 27T 


(5.23) 


is  known  as  the  upwash  or  the  lift  interference  factor. 

In  order  to  include  terms  containing  higher  powers  of  the  ratio  c/h  we  would  have  to  carry 
the  higher  order  terms  in  the  Taylor  expansions  of  functions  H0  and  H7  or  , if  using  the  far  field  approx- 
imation, the  higher  terms  of  the  Laurent  expansion  of  F- . It  can  be  shown  that  in  this  way  the  influence 
of  the  pitching  moment,  the  cross-sectional  area  of  the  airfoil,  and  the  higher  moments  of  the  pressure 
and  the  thickness  distributions  can  be  taken  into  account.  However,  if  the  distributions  of  a and  y are 
known,  as  in  the  computations  by  the  present  method,  one  can  directly  use  Equation  (5.19). 

The  corrected  lift  and  drag  coefficients  C 'L  and  Cp  are  obtained  from  the  wind  tunnel  values 
CL  and  CD  by  resolving  the  aerodynamic  force  into  the  directions  parallel  and  normal  to  the  corrected 
direction  of  the  wind  tunnel  stream 


C 


L 


CL  cos  A a 


CD  sinAa 


(5.24) 


Cjj  = CD  cos  Aa  + CL 


sin  Aa 


(5.25) 


- 21  - 


It  is  perhaps  worth  mentioning  that  the  drag  coefficient  C|} , obtained  as  the  nondimensional  component 
of  the  aerodynamic  force  parallel  to  the  wind  tunnel  axis,  is  as  a result  of  wall  interference  is  in  general 
non-zero  even  if  Q = 0. 


The  angle  of  attack  correction  has  no  direct  effect  on  the  pitching  moment,  that  is 


Cm  = 


(5.26) 


The  usefulness  of  the  angle  of  attack  correction  in  practical  applications  is  examined  in  Section  7.0. 
6.0  VORTEX  PANEL  METHOD  FOR  THE  WALL  INTERFERENCE  CASE 


Using  the  results  of  Section  5.0,  it  is  a simple  matter  to  modify  the  vortex  panel  algorithm 
developed  in  Section  4.0  to  include  the  wall  interference  effect. 


In  the  presence  of  wind  tunnel  walls,  Equation  (4.2)  is  replaced  by 


fj+l  f(f)  fj+l  r -i 

/ T~-z  ^ + / L0(f)H°(z>f)  + T(f)H7(z,f)]  Idfl  <61> 

fi  ^ fj 


This  follows  from  Equations  (5.3)  to  (5.5).  In  order  to  evaluate  the  additional  integral,  we  note  that 
the  functions  H0  and  H7,  given  by  Equations  (5.8)  and  (5.9),  are  analytic  in  the  entire  wind  tunnel 
interior  in  both  arguments  z and  f,  and  can  thus  be  expanded  in  the  following  Taylor  expansions  about 
the  panel  midpoint  Zj : 


AjF-(z)  = 


27Ti 


H(z,r)  = H(z,zj)  + 


dH(z,n 

dr 


(r-  Zj)  + ... 


Assuming  that  the  length  ir,+ , - fj|  of  the  panel  is  small  compared  to  its  distance  from  the  wall 
'.!?  lm  zi  ’ we  can  restr*ct  ourselves  only  to  the  constant  term  in  the  above  expansion.  This 


fj+i  1 

/ a(r)H0(z,r)idri  = - <0j  + 0j+1  )H0(z,Zj)irj+l  - ^ 

tj 


(6.2) 


fj+l 

f 7(r)H7(z,r)idri 


1 

2 (yj  + 7j+i  )H7(z,zj)irj+1  - jji 


(6.3) 
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Introducing  complex  coefficients 


Sk.i  = 


- H0(zk,zj)  |fj+1  - , j = m(E-l)  + 1 


-H0(zk)Zj)|fj+1  - fjl  + 


H0(zk,Zj_  ,)!$;-  fj_i  ! 


, j 4 m(C-l)  + 1,  m(8) 


Gk,j  “ 


-H0(zk,Zj_,)|?j-fj_||  , j = m(8) 


- ~ H7(Zkzj)  l?j+i  " ?jl  . j = m(8-l)  + l 


H7  (zk  ,Zj ) If j+  ] - fj  I - 


H^Z^Zj.,  ) Ifj  fj_, 


, j * m(C-l)  + 1,  m(8) 


- H7(zk,Zj_j ) |fj  - fj_,  | , j = m(8) 


and  substituting  Equations  (6.1)  to  (6.3)  in  Equation  (4.1),  we  obtain 


m (n)  m (n) 

F-(zk)  = 2 (Ck>j  +Skij)aj  + i 2 (Ck (j  + Gk j ) 7j 

j ■ l j - l 


which  replaces  Equation  (4.9).  As  before,  the  coefficients  Ck  -}  are  given  by  Equations  (4.10)  to  (4.12). 

Since  the  oncoming  wind  tunnel  stream  is  at  angle  a = 0 with  the  x-axis,  the  airfoil  boundary 
condition.  Equation  (4.13),  takes  the  form 


Re  | [1  + F"  (zk ))  e^k  | ^ (ok  +ok  + 1) 


Substituting  for  F~  (zk  ) from  Equation  (6.6)  we  again  obtain  a set  of  linear  algebraic  equations  described 
by  Equation  (4.14)  except  that  now 

AgMj  " to[e‘rk(qlj*Gk,,)]  . k*m(«-l),m(«)  (6.8) 
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and 


} 

l 


-)k+ 1 Re 


m (n)  -j  | 

(CM  + Sk.j)oj]J 


1 

2 <ak  + °k  + j ).  k t m(fi-l),  m(C)  (6.9) 


The  elements  of  the  remaining  rows  follow  from  the  trailing  edge  conditions  and  are  again  given  bv 
Equations  (4.17)  to  (4.20).  J 

The  overdetermined  system  of  linear  algebraic  equations  is  solved  as  described  in  Section  4.0 
and  the  values  of  the  airfoil  pressure  coefficients  are  obtained  from  Equation  (4.23).  The  pressure 
coefficients  along  the  wind  tunnel  walls  can  be  evaluated  from  Equations  (4.1)  and  (6.1)  to  (6.3) 
using  the  Bernoulli  theorem:  ’ 


Cp(z)  = 1 - |1  + F-(z)|2 

where 


(6.10) 


Im  (z)  = ± h/2 


In  comparison  with  experimental  wall  pressure  coefficients,  they  can  be  used  to  provide  a posteriori 
check  on  the  values  Py  and  PL  which  were  used  in  the  above  computation.  A practical  implementation 
of  this  technique  has  been  reported  in  Reference  4. 


7.0  NUMERICAL  EXAMPLES 


In  order  to  demonstrate  the  accuracy  and  efficiency  of  the  present  vortex  panel  method, 
example  computations  are  performed  for  some  theoretical  profiles,  whose  exact  potential  flow 
solutions  are  known.  The  exact  solutions  assume  tangent  flow  conditions  all  over  the  airfoil,  so  that 
stagnation  points  are  obtained  at  finite-angle  trailing  edge. 


First,  a circulatory  flow  past  the  unit  circle  is  evaluated.  The  corner  points  are  selected  at 


and  a 10  degree  angle  of  attack  is  assumed,  a = it/ 18.  Table  1 compares  the  exact  pressure  coefficients 
with  the  L.S.  (least  squares),  L] , and  L„  solutions.  It  is  observed  that  the  agreement  with  the  exact 
solution  is  excellent,  and  that  all  three  methods  give  practically  identical  results.  From  the  technical 
point  of  view  it  does  not  seem  to  make  much  difference  which  type  of  the  residual  norm  is  minimized. 
As  far  as  the  cost  is  concerned,  the  execution  time  for  the  single  precision  vortex  panel  method  on  the 
IBM  370/3032  computer  was  0.5  sec.  for  the  L.S.  solution,  0.7  sec.  for  the  L,  solution,  and  0 6 sec 
for  the  L„  solution. 

Table  2 presents  the  results  of  similar  computations  on  the  cambered  Karman-Trefftz  airfoil, 
whose  co-ordinates  are  obtained  using  the  mapping  function 


p = 179/90,  s = 0.92 
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from  the  unit  circle  co-ordinates 


fk  = s + - l]  e~‘£,  e = tt/36,  k = 1,  2, . . 37 


By  the  Schwarz-Christoffel  theorem,  the  theoretical  trailing  edge  angle  is  0 = tr(2  - p)  = tt/90,  that  is 
two  degrees.  The  flow  is  assumed  to  be  at  zero  incidence  with  respect  to  the  real  axis,  a = 0.  Compared 
to  the  circle,  the  agreement  between  the  exact  and  computed  pressure  distributions  on  the  Karman- 
Trefftz  airfoil  is  not  as  good,  notably  in  the  vicinity  of  the  sharp  trailing  edge.  Substantial  improvement 
in  accuracy  can  be  achieved  by  doubling  the  number  of  panels  as  demonstrated  in  Table  3,  but  of  course 
at  the  increased  cost  of  computation.  The  L.S.,  Lj , and  L„  solutions  took  2.5,  4.0,  and  2.9  sec.,  re- 
spectively. However,  from  Figure  5,  showing  the  exact  and  the  L.S.  solutions  in  the  unit  chord  scale, 
it  may  be  seen  that  with  72  panels  the  agreement  is  more  than  satisfactory. 

As  the  third  test  example  is  chosen  the  two-component  airfoil  by  Williams,  Reference  22. 

The  contour  points  of  the  “Configuration  A”,  tabulated  in  Reference  22,  are  directly  taken  as  the 
corner  points  for  the  present  vortex  panel  method.  Table  3 compares  the  pressure  coefficients  by 
Williams  (CP  EXACT)  with  the  computed  ones  by  the  present  method,  using  the  least  squares  mode 
(CP  L.S.).  From  the  pressure  curves  in  Figure  6 it  is  seen  that  the  vortex  panel  method  predicts  cor- 
rectly the  suction  peaks  on  both  airfoil  components  and  that  the  overall  agreement  is  excellent  except 
in  the  trailing  edge  region  of  the  main  airfoil.  In  conjunction  with  vortex  panel  methods,  a similar 
trailing  edge  discrepancy  on  the  Williams’  airfoil  was  also  observed  by  Choo,  Reference  23,  but  appar- 
ently not  by  Weibust,  Reference  24,  who  did  not  force  the  flow  through  stagnation  points  at  the 
trailing  edges,  assuming  that  they  behaved  like  cusps.  However,  since  the  Williams’  pressure  coefficients 
were  verified  to  within  a three  digit  accuracy  by  Ives,  Reference  25,  using  a modernized  Theodorsen- 
Garrick  conformal  mapping  technique  for  the  analysis  of  two-component  airfoils,  it  is  safe  to  assume 
that  the  observed  trailing  edge  discrepancy  is  the  liability  of  the  present  vortex  panel  method.  Using 
this  and  the  preceding  example,  one  may  conjecture  that  it  is  difficult  to  correctly  reproduce  large 
pressure  gradients  near  trailing  edges  unless  the  number  of  panels  is  increased  significantly.  The 
mapping  of  the  trailing  edge  onto  a circular  segment  is  by  its  nature  much  better  suited  for  the  task. 
Nevertheless,  the  main  attractiveness  of  the  panel  method  is  its  simplicity  and  direct  applicability  to 
airfoils  having  a number  of  components  larger  than  two,  where  the  method  of  Reference  25  cannot  be 
used,  and  to  complex  geometric  configurations,  such  as  multi-component  airfoils  between  wind  tunnel 
walls.  Apparently,  one  could  exploit  the  advantages  of  both  techniques  by  removing  sharp  trailing 
edges  (or  edges  having  a small  radius  of  curvature)  by  Karman-Trefftz  transformations  and  applying 
the  panel  method  to  the  transformed  configuration.  However,  since  complications  arise  with  regard 
to  the  present  treatment  of  the  wall  interference  effect  and,  in  addition  to  it,  the  conformal  mapping 
techniques  seem  to  fall  somewhat  outside  the  scope  of  the  present  report,  no  example  is  given  here. 

To  obtain  the  overall  comparison  between  the  computed  and  the  theoretical  results, 

Table  5 presents  the  lift  and  drag  force  coefficients,  obtained  by  the  integration  of  pressure  distri- 
butions. Compared  are  the  exact  values  and  the  computations  by  the  source  panel  method  (Douglas- 
Neumann  type)  from  Reference  22  with  the  computations  by  the  present  vortex  panel  method.  The 
integration  of  pressure  coefficients  over  the  co-OTdinate  points  of  Table  4 is  performed  using  the 
method  by  Nonweiler,  Reference  26.  The  force  coefficients  are  normalized  by  the  chord  length 
c = 1.33.  It  is  observed  that  the  present  method  predicts  the  lift  forces  with  a better  accuracy  than 
the  source  panel  method,  but  more  importantly,  the  resultant  drag  force  is  only  0.0006,  which  is 
acceptably  close  to  the  zero  value  according  to  the  d’Alembert  paradox.  This  result  indicates  that 
the  present  method  is  accurate  enough  to  be  used  in  conjunction  with  viscous  calculations  for  the 
evaluation  of  drag. 

The  column  CPW  L.S.  in  Table  4 gives  the  pressure  coefficients  calculated  by  the  present 
method  (in  the  least  squares  mode)  for  the  following  wall  interference  case:  c/h  = 0.25,  = 1.2, 

PL  = 0.8.  As  before,  the  airfoil  is  at  zero  incidence  with  respect  to  the  direction  of  stream  at  infinity 
(the  wind  tunnel  axis).  Due  to  the  wall  effect,  the  lift  force  decreases  from  its  free  air  value  of 
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2.7887  to  2.2854  and  the  drag  force  increases  from  0.0006  to  0.1598,  cf.  Tables  5 and  6.  In  order 
to  clarify  the  major  trends  of  these  changes,  Equation  (5.21)  is  used  to  evaluate  the  angle  of  attack 
correction  A a = - 3.958  deg.  and  the  lift  and  drag  coefficients  are  corrected  according  to  Equations 
(5.24)  and  (5.25).  From  Table  6 it  is  seen  that  the  correction  value  2.2910  does  not  greatly  differ 
from  the  uncorrected  value  2.2854,  but  the  effect  on  the  drag  coefficient  is  substantial:  the  corrected 
value  0.0016  is  again  very  close  to  zero.  For  comparison,  in  column  CPF  L.S.  on  Table  4 printed  are 
pressure  coefficients  of  the  free  air  case  at  the  angle  of  attack  a = -3.958  deg.  The  values  of  the  lift 
and  the  drag  coefficients,  obtained  by  the  integration  of  pressure  coefficients,  are  given  in  Table  6. 

The  departure  of  the  free  air  value  2.3468  of  the  lift  coefficient  from  the  corrected  wind  tunnel  value 
2.2910  can  be  verified  on  the  difference  of  pressure  curves  in  Figure  7.  The  corresponding  drag  loops 
are  displayed  in  Figure  8.  Here  the  pressure  coefficients  on  the  main  airfoil  and  the  flap  plotted 
against  the  co-ordinate 


yn  = (y  - y, ) cos  a - (x  - xt ) sin  a 


where  xt , y,  is  the  respective  trailing  edge.  It  is  observed  that  in  free  air  the  negative  contribution  of 
the  main  airfoil  is  counterbalanced  by  the  positive  contribution  of  the  flap.  In  the  wind  tunnel  case 
the  flap  contributes  considerably  more  than  the  main  airfoil,  but  the  balance  of  drag  contributions  is 
regained  by  applying  the  angle  of  attack  correction  to  the  co-ordinate  yn . However,  the  remaining 
differences  between  the  pressures  in  the  free  air  case  and  the  corrected  wind  tunnel  case  which  are 
observed  on  upper  surfaces  of  both  the  main  airfoil  and  the  flap  are  again  of  the  same  magnitude  as 
those  seen  in  Figure  7.  This  indicates  that  the  wind  tunnel  case  corrected  for  Aa  and  the  free  air  case 
are  not  exactly  the  same*.  The  difference  between  the  two  flow  patterns  is  also  apparent  on  the 
pressure  distributions  computed  at  the  location  of  the  walls,  Figure  9.  Although  the  airfoil  in  the  wind 
tunnel  at  a = 0 deg.  and  in  free  air  at  a = - 3.958  deg.  produces  nearly  the  same  amount  of  lift,  the  wall 
pressure  signatures  differ  considerably. 

A more  direct  approach  is  to  avoid  the  concept  of  the  angle  of  attack  correction  altogether 
and  interpret  the  corrections  to  pressures,  force  and  pitching  moment  coefficients,  and  other  aero- 
dynamic quantities  of  interest  as  the  differences  between  their  free  air  and  wind  tunnel  values, 
computed  at  the  same  angle  of  attack.  For  example,  as  the  corrections  to  pressure  coefficients  on  the 
Williams’  airfoil  we  can  use  the  differences  between  the  values  CP  L.S.  and  CPW  L.S.  in  Table  4. 
Another  example  is  given  in  Figure  10  where  the  lift  coefficient  CL , the  drag  coefficient  Co . ar>d  the 
quarter-chord  pitching  moment  coefficient  CM  are  evaluated  on  the  Williams’  airfoil  in  free  air  and  the 
wind  tunnel,  at  angles  of  attack  - 6,  - 4,  - 2,  0,  and  2 deg.  Using  the  interpolation  curves,  corrections 
to  CL , CD , and  CM  can  be  found  on  the  whole  interval  - 6 deg.  < a < 2 deg.  Naturally,  this  approach 
can  be  of  practical  value  for  correcting  the  wind  tunnel  data  only  after  the  viscous  effects  on  the  air- 
foil have  been  incorporated  and  the  selected  values  of  porosity  parameters  checked  by  comparing  the 
experimental  and  computed  wall  pressure  distributions. 

For  comparison,  in  Figure  10  are  as  dash-dotted  lines  also  drawn  the  wind  tunnel  values  of 
CL , CD , and  CM  at  a = - 4,  - 2,  0,  and  2 d corrected  according  to  Equations  (5.24)  to  (5.26)  for  a 
variable  (negative)  Aa.  It  is  seen,  for  example,  that  the  line  of  corrected  CL  at  a = 0 deg.  intersects 
the  line  of  free  air  CL  near  a - - 4.4  deg.,  but  the  corrected  wind  tunnel  CD  at  a = 0 deg.  agrees  with 
the  free  air  CD  near  a = - 4.0  deg.  — a difference  of  0.4  deg.  Apparently,  with  respect  to  the  force 
coefficients  the  angle  of  attack  correction  is  a tool  suffering  from  a great  deal  of  astigmatism  in  this 
particular  case:  we  can  either  focus  on  the  free  air  CD  or  on  the  free  air  CL , but  not  at  both  at  once. 
As  far  as  the  pitching  moment  is  concerned,  the  situation  is  even  worse:  we  are  completely  out  of 
“focus”  in  the  given  interval  of  angles  of  attack.  (In  absolute  terms,  however,  the  difference  between 
the  two  sets  of  CM  values  is  quite  small.) 


From  the  engineering  point  of  view,  however,  these  differences  may  be  considered  insignificant. 
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It  should  be  realized  that  the  studied  wall  interference  example  is  somewhat  arbitrary, 
although  by  no  means  unrealistic.  However,  it  serves  to  illustrate  that  the  simple  angle  of  attacx  cor- 
rection procedure  may  have  certain  limitations,  especially  when  applied  to  high  lift  systems. 

As  pointed  out  above,  since  the  present  method  can  be  used  to  compute  both  the  free  air 
and  the  wind  tunnel  cases  at  the  same  angle  of  attack,  the  difference  of  these  sets  of  data  can  directly 
be  used  as  correction  to  the  measured  wind  tunnel  data.  This  process  would  eliminate  the  need  for 
angle  of  attack  correction.  Furthermore,  in  conjunction  with  the  modelling  of  displacement  effects 
the  present  method  also  contains  the  mechanism  by  which  the  wall  effects  on  the  development  of 
boundary  layers  can  be  accounted  for  — something  that  cannot  be  done  with  conventional  wall  cor- 
rection techniques. 

On  the  IBM  370/3032  computer  the  execution  time  required  to  compute  the  122  panel 
Williams’  airfoil  was  9.6  sec.  in  the  free  air  case  and  17.4  sec.  in  the  wind  tunnel  case.  The  increase 
of  the  computation  time  in  the  wind  tunnel  case  is  due  to  the  fact  that  the  elements  of  the  matrix 
carry  additional  terms  which  require  the  computation  of  Green’s  function  values.  In  the  case  that  the 
inverse  of  the  matrix  has  been  generated,  the  repeated  computation  for  a different  angle  of  attack  takes 
only  about  0.2  sec.  in  the  free  air  case.  For  the  wind  tunnel  case,  unfortunately,  any  change  of  angle 
of  attack  gives  rise  to  a completely  new  airfoil-wind  tunnel  configuration  so  that  a new  matrix  has  to 
be  assembled  each  time.  However,  for  repeated  computations  where  only  the  normal  component  of 
velocity  is  affected  (iterations  on  the  boundary  layer  displacement  effect),  the  matrix  of  the  system 
does  not  change  so  that  the  subsequent  use  of  its  inverse  brings  about  the  same  time  economy  as  in 
the  free  air  case. 


8.0  CONCLUSIONS 

The  developed  vortex  panel  algorithm  for  the  computation  of  potential  flow  past  two- 
dimensional  airfoils  is  very  simple  and  efficient,  and  at  the  same  time  capable  of  providing  sufficient 
accuracy  for  most  practical  applications.  For  a further  improvement  in  the  reproduction  of  large 
pressure  gradients  near  sharp  trailing  edges,  or  edges  having  a small  radius  of  curvature,  the  use  of  the 
present  method  in  conjunction  with  conformal  mapping  techniques  is  recommended.  However,  even 
in  its  present  form  it  yields  ample  accuracy  to  permit  analysis  of  such  effects  as  the  inviscid  drag  in- 
duced by  porous  wind  tunnel  walls.  It  is  believed  that  the  present  method  will  help  to  accomplish 
the  main  objective  — the  prediction  of  the  aerodynamic  performance  of  high  lift  systems  in  free  air 
and  in  the  wind  tunnel,  with  the  consideration  of  viscous  effects. 
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THIS  PACES  IS  BEST  QUALITY  F1UCT1CA&IS 
FROM  OOPY  nmiSHBD  TO  DDG  — 


TABLE  1 


EXACT  AND  COMPUTED  PRESSURE  COEFFICIENTS  FOR  THE  UNIT  CIRCLE, 

a = 10  DEG. 


1 . 00000 
0. 98481 
0.93969 
0.86603 
0.76604 
0.64279 
0.50000 
0.34202 
0 . 17365 
0.00000 
-0.  17365 
-0.34202 
-0 . 50000 
-0.64279 
-0.76604 
-0.86602 
•0.93969 
-0.98481 
-1 .00000 
-0.98481  - 

-0 ! 7^605  - 

-0.64279  - 

-0.50000  - 
-0.34202  - 

-0.17365  - 

-0.00000  - 
0.17365  - 

0. 34  202  - 

0.50000  - 

0.64279  - 

0.76604  - 

0.86602  - 
0.93969  - 

0.98481  - 

1.00000  - 


0.00000 
0.  17365 
0.34202 
0 . 50000 
0.64279 
0 . 76604 
0 . 86  60  2 

1 .00000 
0.98481 
0.93969 
0.86603 
0.76605 
0 . 64279 
0.50000 
0. 34202 

0 . 17365 
0 . 00000 
0.  17365 
0.34202 
0.50000 
0.64279 
0.76604 
0.86602 
0.93969 
0.98481 

1 . 00000 
0.98481 

0.93969 

0.86603 
0. 76605 
0.64279 
0.50000 
0. 34202 
0.17365 
0.00000 


CP 

EXACT 

1 . 00000 
0.87939 
0.51754 
0.06365 
0.81521 
1 . 66627 

2 .53208 
3.32368 

4.50980 

4 . 36808 
3.9581  1 
3.32368 

2.53209 
1 .66628 
0.81521 
0 . 06366 
0.51754 
0.87938 

1 .00000 
0.88661 
0.57398 
0. 11964 

-0.40373 
-0.91754 
-1 . 34729 
-1  . 63192 
-1.73143 
-I.63192 
-1 . 34730 
-0.91755 
-0.40374 
0.119b3 
0.57397 
0.88660 

1 .00000 


1 . 00000 
0.87936 
0.51752 
-0.06386 
-0.81538 
-1  . 66657 
-2.53231 
-3.32403 

-4.51067 

-4.36821 

-3.95868 

-3.32448 

-2.53232 
-1  . 66643 
-0.81537 
-0.06383 
0.51751 
0.87937 

1 . 00000 
0.88659 
0.57389 


0.11956 

■0.40382 


-0.91770 
■1 . 34769 
-1 .63212 
-1  . 73187 

■1.63198 
■1  . 34767 
■O.9179O 
-0.40383 

0.11958 

8:1  ml 

1 .00000 


1 . 00000 
0.87937 
0.51754 

0 . 06384 
0.81535 
1 . 66652 
2.53226 
3-32397 

3:!18K 

4.51060 
4. 36812 
3.95882 
3.32443 


■0.81534 

-O.O6381 

0.51752 
0.87938 

1 . 00000 
0.88658 
0.57367 
0.  1 1954 

O.40387 

0.91776 
1.34777 
1 . 63219 

1.73194 
1 . 63206 
1 . 34774 
0.91795 
O.40366 
0. 11955 
0.57392 
0.88657 

1 . 00000 


1 . 00000 
0.87936 
0.51754 

O.O6385 
0.81536 
1 . 66656 
2.53227 
3. 32401 

2:11111 
4 .51062 
4.36812 

1:!1S21 

?:HK1 


zUllIl 

-0.81534  -0.81537 

-O.O6381  -O.O6363 

0.51752  0,51752 

0.87938  0.67938 

1.00000  1.00000 
0.88658  0.88658 

0.57367  0.57368 

0.11954  0.11955 

-O.40387  -O.40385 

-0.91776  -O.91774 


-1 . 34774 
-1 .63216 
-1 .73194 
-1  .63205 
-1  . 34769 
-O.9179I 
-0.40366 

8:ip 

0.86658 

1 .00000 
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TABLE  2 

EXACT  AND  COMPUTED  PRESSURE  COEFFICIENTS  FOR  THE 
KARMAN-TREFFTZ  AIRFOIL,  a = 0 DEG. 


1 . 82978 
1.79734 
1 . 7047  1 
1 . 55907 
1 . 36828 
1 . 1 4044 
0.88371 
0. 60604 
0.31521 
0.01871 
-0.27623 
-0.56267 
-0.83402 
-1  . 08406 
-1  . 30704 
-1 . 49768 
-1  .65128 
- 1 . 76372 
-1.83155 
-1 .85194 
-1 . 82288 

-1.43675 


•1.21347 
-0.94883 
-0. 64893 
-0.32191 
0.02193 
0.37027 
0.70921 
1 .02417 
1.30106 
1.52761 
1.69446 
1 . 79589 
1 . 62976 


0.00000 

0.00669 

0.02691 

0.05954 

0.10l8l 

0. 14990 
0. 19943 

0 . 24603 
0 . 28572 
0.31519 
0.33205 
0.33491 
0.32342 
0.2982  ' 


28 

26110 
21436 
16120 
10521 
05023 
00000 
04228 
,.07412 
0.09368 
0.10012 
0.09383 
0.07658 
0.05143 
0.02247 
0.00571 
0.02872 
0.04323 
0.04766 
0.04261 
0.03082 
0.01657 
0.00478 
0.00000 


CP 

c,X  ACT 

1 . 00000 
0. 12589 
0.04118 
0.04865 
0 . 1 4420 
24355 


-0 


71139  -0.70899 

64022  -0.63762 


67732 
72358 
74724 
74464 
1 1 
0 

51422 
28281 
0.23061 

1 . 00000 
0.  17135 
0. 16487 
0.21069 
°.  158J5 
0. 08660 

0.01164 
0.05900 
0.12178 
0 . 17470 
0.21665 
0.24704 
0.26563 
0.27248 
0.26796 
0.25306 
0.23100 

1 . 00000 


CP 

L . S. 

1 . 00000 
0. 16534 
0 .04757 
0.03911 
0. 13021 

0.22932 
0.33315 
0.43427 
O.52736 
0 . 60853 
0.67427 
0.72100 
0.74479 
0.74244 


-0.51252 
-0.28675 
0.20538 
0 .99996 
0.  1 3566 
0.  19570 
0.21314 

0.  16039 

0.08814 

0.01305 

0.05782 

0. 12093 

0 . 17448 
0.21747 
0.2500' 
0.27201 


0.28091 
0.27415 
0.25588 
0.24992 
1 . 00000 


CP 

LI 

1 . 00000 
0. 16561 
0.04806 
0.0385  1 
1 295  1 
22852 
3226 

5 1112 

609  1 8 
67500 
0.72147 
0.74507 
0.74255 
O.70898 
0.63749 
0.51230 
0.28643 
0.20578 
99995 
1 3528 
19598 
21338 
, . 16065 
0.08846 
0.01348 
0.05721 
12010 
17195 
21975 
25089 
27282 
28154 
. .27470 
0.25633 
0.25019 

1 . 00000 


CP 
L 00 

1 . 00000 
0. 16357 
0.04645 
0 .03984 

0. 13070 

0.22964 

8:31318 

8:3§I3? 

0. 67420 
0.72090 
0.74467 
0.74230 
0.70885 
0.63750 
0.51244 
28671 
20536 
99996 
1 3567 
19559 
21300 
16024 
08799 
01292 
05793 
0.12100 
0.  17451 

0.21744 
0.24993 
O.27168 
0.28057 
0.27359 
0.25494 
0.24826 

1 . 00000 


THIS  PAGE  IS  BBS!  QUALITY  PRACtXCABfi! 
FROM  COPY  FURMISHED  TO  DDC 


TABLE  3 


EXACT  AND  COMPUTED  PRESSURE  COEFFICIENTS  FOR  THE 
KARMAN-TREFFTZ  AIRFOIL,  a = 0 DEG. 

(a)  Upper  Surface 


1 .82978 
1.82150 
1 . 797  34 
1.75812 
1 .7047  1 
1 .63804 
1.55907 
1 .46882 
1 . 36828 
1 . 25848 
1.14044 
1.01518 
0.6837 1 
0 . 74700 
0.60604 
0.46180 
0.31521 
0.  1672  1 
0.01871 
-0. 12940 
-0 . 27623 
-0.42093 
-0.56267 
-0.70063 
-0.63402 
-0.96208 
-1 . 08406 
-1.19928 
-1  . 30704 
-1 .4067  1 
-1.49768 

-1.71288 

-1.76372 

-1.80341 


•1 . 83155 
•1  . 84782 
•1 .85194 


0.00000 

0.00166 

O.OO669 

0.01514 
0.02691 
0.04181 
0.05954 
0 . 07969 
0.10181 
0 . 12539 
0 . 14990 
0 . 17476 
0. 19943 
0.22336 
0 . 24603 
0.26696 
0.28572 
0.30190 
0.31519 
0.32530 
0.33205 

8:ll*|T 

0.33094 
0.32342 
0.31248 
0.29828 
0.28106 
0.261  10 
0.23875 
0.21436 
0. 18837 
0. 16120 

0.07735 

0.05023 

0.02430 

0.00000 


CP 

EXACT 

1 .00000 
0 . 16979 
0.  12589 
0.08391 
0.04118 
0.00296 
C. 04865 
0.09579 
0.14419 

0. 19356 
0.24355 
0.29377 
0.34377 
0.39311 

0.44130 

0.48782 

0.53216 

0.57379 

0.61220 

0.64688 

0.67732 

O.70305 

0.72358 

0.73846 

0.74724 

0.74946 

0.74464 


0.64022 

0.58591 

0.51423 

0.41780 

0.28281 

0.08253 


0.23061 
0.68687 
1 . 00000 


1 . 00000 
0.  19699 
0. 13090 

0 . 08692 
0.04565 
0.00296 
-0.04218 
-0.08955 

-0. 13889 
-°.  18933 
-0.24049 
-0.29107 
-0.34155 
-0.39115 
-0.43974 
-0 .48631 

-0.53100 
-0.57255 
-0.61128 
-0. 64553 
-0.67673 
-0.70198 
-0 .72278 
-0.73770 
-0. 74620 
-0.74878 
-0.74366 
-0.73155 
-0.71018 
-0 . 68068 
-0.63912 
-0.58502 
-0.51349 
-0.41747 
-0.28325 
-0.08539 
0.22338 
0.68002 
1 . 00000 


1 . 00000 
0. 19700 
0 . 1 3090 
0.08694 
0.04569 
0.00302 
0.0421  1 
0.08942 

0. 13878 

0.  18924 
0.24039 
O.29098 

0.34144 

0. 39193 

0.43961 

0.46618 

0.53088 

0.57252 

0.61136 


1 . 00000 
0.  1 9649 
0.13062 
0 . 08674 
0.04552 
0.00286 
-0.04208 

-0 . 08942 
-O.I3679 
-0 . 1 8925 

-0.24043 

-0.29102 

-0.34152 

-0.39116 

-O.43975 


-0.64565 
-0.67684 
-0.70207 
-0.72286 
-O.73776 
-0.74625 
-0.74881 
-0.74368 
-0.73157 
-0.71018 
-C  . 68067 
-O.639IO 

-0.58499 
-0.51346 
-0.41742 
-0. 28320 
-0.08534 
0.22344 
0.68006 
1.00000 


0.48628 
0.53093 
0 . 57249 
0.61123 
0. 64552 
0.67673 
0.70197 
0.72278 
0.73768 
0.74618 
0.74876 
0.74363 
0.73152 

0. 7I015 

0.68064 

0.63908 
0 . 56497 

0.51344 
0.41741 
0.28320 
0.08533 
0.22357 
0.67990 
1 . 00000 
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TABLE  3 

EXACT  AND  COMPUTED  PRESSURE  COEFFICIENTS  FOR  THE 
KARMAN-TREFFTZ  AIRFOIL,  a = 0 DEG. 

(b)  Lower  Surface 


■ 1 .8519'* 
-1  . 85307 
- 1 . 82288 
-1 . 78952 
-1  - 74358 
-1  . 68515 

■ 1 .611438 
-1 .53158 
-1  . 43675 
■1 . 33059 
■1.21357 
■1 .08596 
■0. 941883 
•0.80282 
■0. 64489  3 
■0. 418822 
•0.32191 
■0.  15137 

0.02193 
0 . 19638 
0.37027 
0.544  1 82 

o. 7092 1 

0.87060 
1 .0241  1 7 
1 . 16820 
1.30107 
1.52129 
1.52761 
1 . 61894 
1 . 69445  6 
1.75357 
1 . 79589 
1 .82129 
1 . 82976 


0 . 00000 
0.0223 1 
0. 041228 
0.05964 
0.07512 
0.0855 1 
0.09368 
O.O9854 
0.  10012 
0.09848 
0.09383 
0.08641 
0.07658 
0.06476 
0.05143 
0.03715 
0.02247 
0.00800 
0.00571 
0.0l6l0 
0.02872 
0.03719 
0.05323 
0.C4672 
0.04766 
0.04620 
0.04261 
0.03732 
0.03082 
0.02370 
0.01657 
0.01007 
0.00476 
0.00127 
0 . 00000 


CP 

EXACT 

1 . 00000 
0.65458 

-8:Mi?3 

-0  . 18486 
-0.21496 
-0. 21069 
-0.  18919 

-0. 15675 

-0.  12372 
-0.08659 
-0 .04890 
-0.01  164 
0.02449 
0.05900 
0.09152 
0.12178 
0.14956 
0 . 17570 
0.  19709 
0.21665 
0.23331 
0.25704 
0.25781 
0.26563 
0.27051 
0 . 27259 
0.27161 
0 . 26796 
0.26170 
0.25306 
0 . 25249 
0.23099 
0.22165 

1 . 00000 


1 . 00000 
0.64499 
0.  15979 
■ 0 . 08656 


0 . 18752 
0.21573 
0.21135 
0.  18925 
0.  15941 

0. 12374 
0.08697 
0.04932 
0. 0 1 197 
0.02404 
0.05880 
0.09110 
0.12169 

0. 14929 
0.  17567 
0.19699 
0.21666 
0.23360 
0.25744 


0.25853 
0.26673 
0.27229 
0.27492 
0.27518 
0.2721 1 
0.26583 
0.25602 
0.25484 
0.23376 
0.23938 
1 .00000 


1 . 00000 
0.64495 
0 . 1 5974 
-0. 08664 
-0. 1 8746 
-0.21576 
-0.21136 
-0.  18927 
-0.15944 
-0  . 1 2378 
-0.08701 
-0.04937 

-0. 01203 

0.02397 
0.05871 
0.09099 
0.12155 
0.  14910 
0. 17398 
0.19767 
0.21687 
0.23375 
0.25757 
0.25865 
0.26683 
0.27238 
0.27502 
0.27528 
0.27218 
0.26589 
0.25607 
0.25487 
0.23377 
0.23939 

1 . 00000 


1 . 00000 
0.64502 
0.  15982 
0.08655 

0. 18739 

0.21569 
0.21  129 
0.  18920 
0.  15929 
0.12383 

0.08699 

0.04933 

0.01198 

0.02403 

0.05879 

0.091 10 
0.12160 
0. 14907 

0.17590 
0 . 19706 
0.21683 
0.233^5 
0.25727 
0.25867 
0.26678 
0.27235 
0 . 27500 
0.27531 
0.27228 
0.26570 
0.25590 

0.25567 

8:11118 
1 . 00000 
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TABLE  4 

EXACT  AND  COMPUTED  PRESSURE  COEFFICIENTS  FOR  THE  WILLIAMS’  AIRFOIL 

(a)  Main  Airfoil 


1 .00000 
0.99753 
0.99043 
0.97850 
0.96122 
0.93815 
0.90929 
0.87501 
0.63597 
0.  -*9290 
o.,t4656 
0.69769 
0.64697 
0.59507 
0.54258 
0.49010 
0.43814 
0.38724 
0.33785 
0.29044 
0.24543 
0.20321 
0.16414 
0.  12657 
0.09681 
0.06914 
0.04562 
0.02707 
0.01311 
0.00409 
0.00017 
0.00143 
0.00783 
0.01927 
0.03560 
0.05663 
08214 
11190 


14563 
‘ j04 


18^ 
22381 
0.26759 
0.31398 


0.36259 

0.41297 

0.46466 


0.00590 
0.00718 
0.01041 
0.01458 
0.01893  - 

0.02323  - 

0.02762  - 

0.03232  - 

0.03743  - 

0.04293  - 

0.04870  - 

0.05453  - 

0.06020  - 
0.06550 
0.07019 
0.07408 
0.07700 
0.07881 
0.07942 
0.07875 
•0.07678 
0.07352 
0.06902 
O.O6335 

0.05663 
0.04898 
0.04056 
0.03155 
0.0221  1 
0.01242 
0.00264 
-0.00714 
-0.01681 
-0.02625 


■0.03530 
■ 0.04378 
■0.05151 
■0.05833 
■0.06409 
■0.06666 
.0.07194 

•0.07386 
0.074] 


38 

0.07350 
0.07124 
_ . . _ 0.06766 

0.51716  -0.06287 

0.56993  -0.  ' 

0.62240 
0.67396 
O.72396 
O.77I69 
0.81639 
0.85728 
0.89358 
0.92461 
0.94998 
n oAgtc 


8:lttu 

0.99417 

0. 99931 

1 . 


■0.0501 
■0.04261 
■0.03451 
■0.02612 
-0.01771 
-0.00965 

■0.00236 

0.00361 

0.00766 

0.00941 

0.00903 

0.00748 

0.00612 

0.00590 


CP 

EXACT 

1.00000 
■0.92119 
•1 .6OI69 
■ 1 . 9537  4 
■1.95727 
1 .8O9OO 
1.65810 
1.55225 
1.49342 
1 • 47265 
1.48106 
1.51204 
1.56109 
1.62525 
1.70272 
1 .79260 
1 . 89477 
2.00992 

htilk 

2.45564 

3.18166 

3.55471 

4.05084 

4.73940 


CP 
L . S . 


-7 


,00000 

0458’7 


1 
0 

-1.28191  -1 

! 

■ 1 


1 .816^8 


881 
7588? 
1.62122 
1 . 52459 
1.47106 
453^; 


46572 
49832 

54760 

ilffl 

109  6 


8.73166 
■8.34989 
■4.20249 
-0.79779 

0. 53705 

0.93296 

0.99969 

0.95724 

O.88385 
0. 80830 

0.74118 

0.68596 

0.64327 

0.61252 

0.59271 

0.56264 

0.58115 

0.58708 

8:gT13S 

0.63896 
0.66437 
0.6920 
0.7205 
0.74766 
0.76914 
0.77689 

O.7547? 

o\ 50^6^ 

0.23687 

-0.04022 

1.00000 


J « 1 J 

1 . 88299 
1 .99672 

2.  12720 
2. 27347 
2 .44087 

I:l?fl8 

3.  15722 
3.52746 
4.01407 
4.68560 
5.66525 
7.04428 
8.63074 

8.23135 

4.06761 

0.72045 

0.56552 

0.94009 

0.99923 

0.95401 

0.67971 

0. 80372 

0.73664 

0.68166 

0.63916 

O.6O673 

0.58900 

0.57916 

0.57623 

O.56434 
0.1 
O.i 
0.63732 
0.66310 
0.69122 
0.72072 
0.74876 
0.77198 

0.76302 

0.76694 
0.69919 
0.56694 
0.53967 
0.92122 
1 . 00000 


CPW 
L.  S. 

, 00000 
,10010 


64761 
70060 
- 1 .57704 

-I.43627 

-1.33714 
-1.27514 
-1.24534 
-1.24185 
-1.25516 
-I.26176 
-1  . 32067 
-1  . 36771 

-1.55112 
-1 . 62717 
-1.70994 
-1 . 60202 
-1 .90580 
-2.62934 
-2.17355 
-2.35716 
-2.59032 
-2.69974 
32f  “ 
161 

.27021 
.3.34768 

.6.85589 

0.62225 

0.98451 

0.96878 

O.8633O 

0.75138 

0.65537 

0.57964 

8:5!! 


0.44590 

0.44408 

0.45141 

0.46665 

0.48763 


59689  0.51385 

61473  0.54409 


0.54409 

0.57794 

0.61382 

0.65101 

0.66673 

0.72417 

0.75422 

0.77172 

0.76243 

0.70342 

0.58296 

0.56155 

0.92522 

1.60000 


CPF 
L* . S • 

1.00000 

c. 08039 

-1.19640 
-1.70556 
-1.75970 
-1.63325 
-1.49139 
-1 . 36770 
-1  . 32463 
-1 .29462 
-1 .29149 
-1 . |050' 


-1 . 37364 
-1.42234 

-I.47939 

-1.54353 

-1.61352 

II  I7784l 
-1.67456 
-I.96304 
-2.11143 
-2.26326 
-2.45386 
-2.69696 

-3!  46p7 
-4.02814 

-4.46058 

-3.51898 

-6.94154 

0.59881 


0.9618 
0.9702 
0.66379 
O.7494- 

0.6505 
0.57168 

9:tm* 


W 0 ~ C.  J I V « T ( I V I 

0.45806  0.44578 


0.43259 

O.43016 

0.43711 

0.45218 

0.47344 

0.50008 

0.53090 

0.56555 

0.60234 

0.64060 

0.67941 

.0.71605 

0.74712 

0.76538 

0.75621 

8:8  US 

8 :3H1Z 

1.00000 


IMJ 


-r— 


: 4\  , 


: v-  ' 
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TABLE  4 

EXACT  AND  COMPUTED  PRESSURE  COEFFICIENTS  FOR  THE  WILLIAMS’  AIRFOIL 

(b)  Flap 


K 


6 
6 

65 

66 

67 

68 

69 

70 

71 

72 

H 

75 

76 

n t 

7 5 

So 

81 

82 

U 

85 

86 

8 
8 

89 

90 

91 

92 


98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
1 1 1 
112 
1 1 
1 1 
11 
116 

1!  I 

123 
121 
122 
12 
12 


1.31389 
1.31303 
1 .31064 

1:1811? 
1 . 29624 
1.28934 
1.28147 
1.27267 
1 .2629& 
1 .25245 
1.24112 
1 . 22907 
1.21637 
1.20310 
1 . 18938 
1 . 17532 
16103 
1 . 14667 

1:11115 

1 . 10447 
1.09119 
1.07852 
1 . 06658 
1.05548 
1.04527 
1.03600 
1.02768 
1.02027 


1 


93  1.01372 

94  1.00797 

15 


1.00295 
0.99864 
0.99508 
0.99242 
0.99087 
0.99073 
0.99226 
0.99572 
1.00134 
1.00930 
1.01972 
1.03263 
1.04799 
1.06565 
1.08535 
1.10676 
1 . 12944 

l:l?i« 

1.24177 

1:i  m 

1.28864 

?:|8I|I 
1.31121 
1 • 3 1 360 
1.31389 


-0 ; 

-0.20221 
-0. 19840 
-0. 19245 
-0.  18476 
-0.17566  - 

-0.16544  - 

-0.15440  - 

-0.14281  - 

-0.13090  - 

-0.11891  - 

-0.10704  - 

-0.09548  - 

-0.08439  - 

-0.07391  - 

-0.06416  - 

-0.05525  - 
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TABLE  5 

EXACT  AND  COMPUTED  FORCE  COEFFICIENTS  FOR  THE 
WILLIAMS’  AIRFOIL  IN  FREE  AIR,  a = 0 DEG. 


Lift  Force  Coefficients 


Main  Airfoil 

Flap 

Total 

Exact 

2.1853 

0.6242 

2.8095 

Source  Panel 

2.1579 

0.6003 

2.7582 

Vortex  Panel 

2.1662 

0.6225 

2.7887 

Drag  Force  Coefficients 


Main  Airfoil 

Flap 

Total 

Exact 

-0.2886 

0.2886 

0.0000 

1 

Source  Panel 

-0.2821 

0.2676 

-0.0145 

Vortex  Panel 

-0.2861 

0.2867 

0.0006 

TABLE  6 


COMPUTED  FORCE  COEFFICIENTS  FOR  THE  WILLIAMS’  AIRFOIL  IN  THE 
WIND  TUNNEL  AND  FREE  AIR,  BY  THE  VORTEX  PANEL  METHOD 


Lift  Force  Coefficients 


Main  Airfoil  Flap  Total 


Wind  Tunnel 
(a  = 0.000° ) 

1.6906 

0.5949 

2.2854 

Free  Air 
(a  = -3.958°) 

1.7200 

0.6268 

2.3468 

Wind  Tunnel 
Corrected  for 

1.6783 

0.6127 

2.2910 

Aa  = - 3.958° 

Drag  Force  Coefficients 


Main  Airfoil 

Flap 

Total 

Wind  Tunnel 
(a  = 0.000° ) 

-0.1187 

0.2785 

0.1598 

Free  Air 
(a  = -3.958°) 

-0.2433 

0.2426 

-0.0007 

Wind  Tunnel 
Corrected  for 

-0.2351 

0.2367 

0.0016 

Aa  = -3.958° 
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FIG.  5:  EXACT  AND  COMPUTED  PRESSURE  DISTRIBUTIONS  ON  THE 
KARMAN-TREFFTZ  AIRFOIL,  a * 0 DEG. 


FREE  AIR  AT  a = -3.958° 
WIND  TUNNEL  AT  a =0.000' 
c/h  = 0.25,  Py  = 1.2,  PL 

(122  PANELS) 


FIG.  7:  PRESSURE  DISTRIBUTIONS  COMPUTED  FOR  THE  WILLIAMS'  AIRFOIL 
IN  FREE  AIR  AND  IN  THE  WIND  TUNNEL 


MAIN  AIRFOIL 


AG  LOOPS  COMPUTED  FOR  THE  WILLIAMS'  AIRFOIL 
IN  FREE  AIR  AND  IN  THE  WIND  TUNNEL 


